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ABSTRACT 



This work addresses the design and analysis of the 
Pitch, Yaw and Roll autopilot for application to the 
Bank-To-Turn (BTT) missiles. At first, the linear 
uncoupled channels were designed and analyzed according to 
the desired requirements. Utilizing the uncoupled 
channels, the linear coupled autopilots were designed, not 
including the inertial, kinematic and aerodynamic 
cross-coupling. Then, the nonlinear CBTT autopilots were 
designed and analyzed, using the linear CBTT (coordinated 
Bank-To-Turn) models, which now have coupled with 
kinematic, inertial and aerodynamic cross-coupling. The 
minimization of the above kinematic and inertial coupling 
and their effects were completed using feedbacks of 
angle-of-attack and rate of angle-of-attack in the Pitch 
autopilot . 
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I. 



INTRODUCTION 



In recent years, the application of Bank-to-turn 
guidance to missiles has generated considerable interest. 
This interest has been motivated by the drag reduction, 
which the Bank-to-turn (BTT) controlled missiles offer, 
over the conventional cruciform roll stabilized, 
skid-to-turn controlled missiles. 

Bank-to-turn steering may provide improved performance 
for missile systems. In addition, the ramjet engine inlet 
configuration is believed to provide a greater range 
capability than other inlet configurations. BTT control is 
required to satisfy the sideslip constraints imposed by 
chin inlets. Despite these facts, there are unanswered 
questions concerning BTT systems stability during homing, 
guidance performance, autopilot and guidance logic design 
and subsystem requirements which must be investigated 
before Bank-to-turn steering can be considered a viable 
method of control for high performance missiles. 

Many missile programs were initiated during the past 
decade to improve their capability via Bank-to-turn 
control. The results have advanced the understanding of 
the different missile subsystems. In the autopilot area, 
many types of autopilots have been found which force the 
missile to roll or bank so that the steering maneuver 
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occurs with the airframe oriented in a specified or 
preferred direction with respect to the incoming airstream. 
This entire class of autopilots may be referred to as 
preferred orientation control autopilots. 

Also, simplified guidance studies which neglect radome 
effects and assume coordinated missile motion have shown 
that CBTT can provide acceptable performance with roll 
rates that are not excessive for autopilot control. These 
studies were made for a medium range area defense mission 
and a long range suppression mission and considered both 
high lift and moderate lift configurations [Ref. 1]. 

The present work reviews the design and analysis of 
the pitch. Yaw and Roll autopilot for application to the 
Bank-to-turn (BTT) missiles developed in [Ref. 2]. At 
first, the linear model of the uncoupled autopilot channels 
will be designed and analysed in the critical areas of 
concern regarding the response of CBTT control for the 
circular and elliptical airframes. Then, the designed 
autopilot channels will be coupled using suitable 
aerodynamic models and a measure of sideslip control is 
obtained by roll command to the linearized CBTT autopilot. 
The relative stability of the autopilot branches and means 
for improving stability will be discussed. In a later 
chapter the non-linear three dimensional model will be 
designed, using a three dimensional aerodynamic model and 
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the same control laws described in the linear studies, 
except for a minor modification to the coordinating branch 
dependence on Angle-of-Attack and also the inclusion of 
Anti-Gravity Bias. Next, the kinematic and inertial 
coupling effects between pitch and yaw dynamics will be 
studied. Transients in maneuver plane acceleration are 
caused in the above cross-coupling effects and are in the 
form of overshoots and undershoots. Coupling transients 
will be reduced by increasing pitch stability. The 
technique of feedbacks of Angle-of-Attack and rate of 
Angle-of-Attack for the pitch autopilot, will be 
considered, in order to decrease kinematic and inertial 
cross-coupling effects. Finally conclusions and 
recommendations for future study will be stated. 
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II. 



LINEAR DESIGN AND ANALYSIS OF UNCOUPLED AUTOPILOTS 



FOR CIRCULAR AND ELLIPTICAL AIRFRAME 



A. GENERAL 

The linear design and analysis started with uncoupled 
autopilot channels, i.e. aerodynamic, kinematic and control 
law cross-coupling between Pitch, Yaw and Roll channels 
were removed. The autopilot design technique was 
classical, using a combination of frequency response and 
root locus techniques [Ref. 2] , providing the range of 
required missile body angular rates, control motions and 
minimizing the influence of aerodynamic variations on 
desired responses. To achieve the desired maneuver plane 
acceleration response for the CBTT autopilot, both Pitch 
and Roll uncoupled autopilots were designed to have a time 
constant of 0.5 sec in order to achieve the desired 
maneuver plane speed of response. The Yaw channel which 
follows the Roll channel motion to produce desired 
coordination (i.e. minimization of sideslip angles), was 
designed to have a more rapid response than the Roll 
channel, having 0.39 sec for circular and 0.36 sec for 
elliptical time constant in the maneuver plane acceleration 
response. The whole design was done for the constant 
flight condition of 60 KFT altitude and mach number 3.95. 
The design discussed here was developed by Arrow [Ref. 2]. 
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In this chapter the airframe configurations, the 
aerodynamic models and control laws for the Pitch, Yaw and 
Roll channel are reviewed. The transfer functions for both 
circular and elliptical configuration are derived. An 
analysis for every channel in regard to the acceleration 
responses, the body angular rates and control surface 
deflection is made. Also, a comparison of airframe 
responses for every channel is made, considering the 
relative stability and required control surfaces. 

A general block diagram, of the channels and the 
command, which are applied to the control laws is shown in 
Figure 2.1. A number of the Figures shown in this report 
have been taken from [Ref. 2] . 

B. AIRFRAME CONFIGURATIONS 

The U.S. National Aeronautics and Space Administration 
maintains a continuing research effort in missile related 
aerodynamics. In this effort monoplanar missile concepts 
have been considered, in part, because of the relatively 
low geometric profile that can be achieved for convenient 
carriage [Ref. 1]. The two Airframe Configurations studied 
in this research were taken from [Ref. 1] and are shown in 
Figure 2.2 and Figure 2.3. The circular cross sectional 
body has a closure ratio Aj^ase-^^max 0*69 and the 
occurs at 68.0% body. The Elliptical has an exact 3:1 
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Bank-to-Turn Autopilot (BTT) 
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Pig. 2.2 Model of Circular Configuration - 1/6 Scale 
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Model of Elliptical Configuration - 1/6 Scale 



elliptical cross section, and its cross section area is the 
same as the area of circular body. The span of both 
airframes is the same. The wing area and span for the 
circular concept were chosen as typical of current 
maneuvering missiles. The wing for the elliptical 
configuration was determined by projecting the elliptical 
body on the circular bodywing planform. The resultant 
exposed wing planform then became the wing for the 
elliptical body. The tails for each design had a dihedral 
of 30° degrees as shown in Figure 2.2. In order for tail 
deflection to be compatible with complex surfaces of the 
afterbody of the elliptical configuration, the tail hinge 
line was skewed such that a 10 degrees deflection measured 
at the body-tail juncture had a resultant 7.04 degrees 
surface deflection, as shown in Figure 2.3 and [Ref. 1], 
Comparison of the elliptical cross section concept with the 
circular cross section concept indicates the following: 

1. About 25 percent more normal force, that is nearly 
independent of angle of attack can be achieved at 
supersonic speeds. 

2. Values of the longitudinal stability parameter 

a 

are more positive, with more pronounced nonlinearities in 
pitching moment occuring at subsonic speeds. 

3. Levels of directional stability are increased and 
are more compatible with levels of longitudinal stability. 
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4. Noticeably more yaw control is available, although 
suitable locations for tails on the body are more limited 
because of the geometry of the elliptical body. 

C. UNCOUPLED LINEAR PITCH CHANNEL 
1 . Aerodynamic Model 

The initial phase in the design of the CBTT 
autopilots involves the design of the uncoupled pitch 
channel. The linear pitch channel dynamic model which is 
used, is shown in Figure 2.5. 

This figure is taken from [Ref. 2]. The above 
model is linearized for stability studies. The following 
three assumptions were made: 

a. Plane X^-Zg is the maneuver plane, shown in 
figure 2.4 

b. Missile is trimmed in pitch, that is My=0, at 
fixed values of q and <5^. 

c. Missile roll rate is constant. Using the 
foregoing assumptions, and the approximation that the 
variation in forward velocity (u), is equal to zero, a 
simplified set of longitudinal equations of motion are 
constructed applicable to short period motion. These 
equations describe the two degree of freedom, short period 
mode longitudinal motion. 

Linearized aerodynamic derivatives for M=3.95 
are provided in Table I. 
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Fig. 2.4 Aerodynamic Sign Convection 
And axis. 
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Fig. 2.5 Uncoupled Pitch Channel Autopilot 



TABLE I 



Linearized Aerodynamic 



Circular 





a 


= 0° 


a= 


10° 


a = 


20° 




- 


.065 


- 


.082 


- 


.111 




- 


.025 


- 


.019 


- 


.003 






0 


- 


. 009 


- 


. 020 


'y 


+ 


.021 


+ 


.022 


+ 


. 028 




— 


.050 


— 


.053 


— 


. 062 


5y 




0 


— 


.016 


— 


.038 


""a 




0 


— 


. 009 


— 


. 022 


C 




0 


+ 


. 018 


+ 


. 044 




+ 


.031 


+ 


. 035 


+ 


.044 


s 

a 


+ 


.15 


+ 


.17 


+ 


. 22 


p 


+ 


.04 


+ 


. 04 


+ 


. 05 


a 


— 


.060 


— 


.065 


— 


.118 






. 080 


' 


.095 - .115 

Reference 
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Derivatives (M = 3.95) 







Elliptical 






a = 


^ 0° 


o = 10° 


a = 


20° 


- 


. 043 


- .054 


- 


.064 


+ 


.024 


+ . 024 




.032 




0 


- .027 


- 


. 040 


+ 


. 016 


+ .015 




.019 


- 


.042 


- .039 


- 


. 045 




0 


- .010 


- 


. 023 




0 


- .006 


- 


. 014 




0 


+ .014 




.032 


+ 


. 023 


+ .023 


+ 


.029 


+ 


. 18 


+ .22 


+ 


. 30 


+ 


.02 


+ .02 


+ 


.025 


+ 


.015 


+ .0137 


- 


. 125 


- 


. 055 


- .055 


- 


.075 


C.G. 


at 0 


.6 
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Pitch Control Law 



The pitch control laws for the circular and 
elliptical airframes are shown in Figure 2.6 as given in 
the [Ref. 2]. Lag-Leads were used to prevent guidance 
noise saturation problems. An integrator was used in the 
acceleration error branch of the circular control law to 
satisfy the guidance requirement of zero steady state 
error. In the elliptical control law, an integrator was 
placed in the rate error path (i.e path after the adder 
(<3+y) / shown in Figure 2.6) for reduction of the 
acceleration response overshoot below 10 percent and a gain 
K=1.1053 was placed in series with the acceleration command 
for zero steady state acceleration error. 

A normal acceleration command n^ equal to 1 gee, 

is applied to the pitch control law which uses measurements 

of missile body pitch angular rate (q) and pitch normal 

acceleration (n_) to determine the required actuator 

command (S ). The actuator is modeled as a first order lag 
^c 

at 188.4 rad/sec and is shown in Figure 2.5. 

3. Transfer Functions of Aerodynamic Models for 

Circular Airframe 

The equations which are represented by the block 
diagram of the aerodynamic model, of the uncoupled Pitch 
channel (Figure 2.5) are: 
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Fig. 2.6 Pitch Control Laws 



S = S'5'dK<a+5T.3C„,4p ] / luu fTaJ/«c] 

^ •’ (II. C. 3-1) 

'yij=-q^[5MQf,-f?+C»,'a]/w [gewj (ii.c. 3 - 2 ) 

dr 57-3[C| + 3a.XTl,/ V] [4e(j] (il.c.3-3) 

Utilizing Laplace Transformation, the above 



equation becomes: 

fr,. ( Cwi^* a + s;.5 Cwjp * cTp j 

jS- - 

lyy 

- 9 ^ (5f-3 a)/w 

5a- 5?.3(q + 



(II. C. 3-4) 

(II.C.3-2) 
(II. C. 3-5) 



Introducing equation II.C.3-2 to II. C. 3-5 the last 
equation yields. 




^ ) 

lyy 



(II.C. 3-6) 



S a - 5^.5 




q ^ (si-i <5p + 
W 




(II.C. 3-7) 

Rearranging above equations, in terms of q and a, 
we obtain: 

Applying the Cramer's rule in the system of 
equation II.C. 3-8 and II.C. 3-9 the transfer functions q/ 
and a/ 6 are obtained. 
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(II. C. 3-10) 



And : 




Where ; 

K = (57.3) (32.2) ^ 

V 2 

S = 3.14159 (ft"^) 

d =2 (ft) 

q = 1650 (Ibs/ft^) 

lyy = 804 (slug-ft^) 

W = 2525 (lbs) 

C„ = -0.06 

a 



= -0*08 
in 




a 



C„ =0.04 
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As referred to in Appendix A and Table I. After 
substituting the aforementioned values into equation 
II. C. 3-10, it becomes: 

^ - - i. 5561 S - 0.^59 



5p 



22.5M5Xld*S^+5.5635X<0?5+i 



( ^Isec) 



(II. C. 3-12) 



Introducing equation II. C. 3-12 into equation 



II. C. 3-2 the last equation yields: 



5T.5K(s>^)-G^pQoc 

W2 



Jk.- 



5- 



3283.29 



lyy • W 






K* q Cwa 



y- 



5T.3 q Cin. 



5T3 q ^ 
W 



W ; Jyy 

( 1 1. C. 3-13) 

Substituting the above mentioned values into 
equation II. C. 3-13 and minor manipulating, it yields: 

^ Z _..iOt OT^X<oV- 19.82 (9«5/racl) 

22. 5^5 3.3633X40 *5 + i (II. C. 3-14) 

The equations II. C. 3-12 and II. C. 3-14 are the 

pitch aerodynamic transfer functions for the pitch angular 

rate q about the Yg and the achieved maneuver acceleration 

in the direction (n„) . 

D Z 

4 . Design Approach and Analysis For Circular Airframe 
A Fixed Flight condition was selected for these 
preliminary performance studies of circular and elliptical 
airframes. The selected flight condition, at 60KFT 
altitude and mach 3.95, provides a sufficiently low dynamic 
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pressure, so that missile maneuvers will result in large 
enough angles-of-attack, to exercise side slip control. 

The above fixed condition was selected, in order to reduce 
the complexities of design of the CBTT autopilots. 



classical, using a combination of frequency response and 
root locus technique [Ref. 2] to achieve practical 
bandwidths (i.e., sufficient high frequency attenuation) 
and therefore provides the range of required missile body 
angular rates and control motions. The analysis of the 
CBTT autopilot channels is based, on the time responses of 
manuever plane accelerartion, body angular rates and tail 
incidence angles. 



written (Appendix D) using the equations which are 
represented by the block diagrams of aerodynamic model and 
linear uncoupled pitch control law, as are shown in figures 
2.4 and 2.6. 

The equations of the uncoupled linear pitch 
aerodynamic model and control law are; 

a. Equations of Control Law 

(1 ) Acceleration Filter Equation 



The uncoupled autopilot design technique was 



For purposes of analysis, a CSMP program was 




Hi 



(II. C. 4-1) 



Rearranging above equation II. C. 4-1 one obtains; 
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(II. C. 4-2) 



X = -150 X+ 

(2 ) Acceleration Compensator Equation 

S 



V= 



(II. C. 4-3) 



Utilizing inverse Laplace Transformation and 
rearranging the above equation, it becomes; 

V + 5H-i.353l(ii2c“X) -■I93.5X^6 ^-Ct11c-X) (ii.c.4-4) 

■jk/ . 

Using state-space representation of a system, 
in which the forcing function involves derivative terms 
[Ref. 5: pp. 675-678] one obtains; 

Yi =Y2 -lass'! ('n?^-x') (II. c. 4-5) 

V 2 r-SYo +G.512 -(n^c-x) (ii.c.4-6) 

Y ■= '/i (II. C. 4-7) 

(3 ) Rate Compensator Equation 




(II. C. 4-8) 



Using Inverse Laplace Transformation, it 
yislds: • , • 

Sp,z -0.<a3'Spc-3M0.5SX^0 (v<-xO-2.230S*(V4-)C() 

^ ^ (II.C.4-8a) 

b. Actuator Equation (Figure 2.5): 

(S'p- ^ (II. C. 4-9) 

— hi 
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or 



<Jp+ (II. C. 4-10) 

c. Equations of Aerodynamic Model 

In order to achieve the plots of the time 

response of q, n , 6 in the linear uncoupled pitch 

z p 

channel, the transfer functions which are represented by 
the aerodynamic model and derived in section II. C. 3, were 
used : 

9 _ ~1»S5U$ — 0.159 (II. c. 4-11) 

25.5i}5Xi0^5^ + 3.3635X^O^S + 1 



or 



^ (25.545X^6^+3.3633X40^5+1) = (- 1.356J -5 - 0.159) • <Jp 

(II.C.4-lla) 

Rearranging and minor manipulating equation 
II.C.4-lla it yields: 

‘<^ + 1M9.1S2X10^C) + M4.5556 ='59.2635 <^p -+.052 -^p (ii.c.4-12) 

Utilizing state representation of a system in 

which the forcing function involves derivative terms, as in 

section II. C. 4, one obtains: 

- X 2 - 59.2635- S'p (II. C. 4-13) 

)(^ Z -i^ii.3556Xi-149.mxm\2+J.?8S<5p (II. C. 4-14) 

q z (II. c. 4-15) 

Also, rearranging and manipulating in the same 

way, equation II. C. 3-12 yields; 

ii - ^2 “■? 0 x ^01 X 10 ^. <Sp (II. C. 4-16) 

r -^4.3356 li - 199.132X10^22 -10433XlO^(fp (ii.c.4-17) 
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+ M.T05 <5p ') (II. C. 4-18) 

Utilizing state-space representation, all the 



aforementioned equations can be modeled in the following 
ninth order system: 



X. 




0 

5 


1 


0 


0 0 


0 


0 


0 


-59.36 








0 ' 


X 2 






-0.15 


0 


0 0 


0 


0 


0 


i.ms 




Xi 




0 


« 

h 




0 


0 


0 


4 0 


0 


0 


0 


-o.m 




2i 




0 


# 




0 


0 


•44.54 


-0.45 0 


0 


0 


0 


-1.043 

- 1 o'l'i.iS' 




22 




0 


• 
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•? ■*' 
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Vi 
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0 


0 0 


1.353 
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0 




Yi 




-i.35 


* 












c> 




















0 


0 


0 


0 0 


6.5T2 


-5 


0 


0 




Ya 




C5f 


^Pc 




2.251 


0.344 


0 


0 0 


-0.46f 


-0.541 


-0.-I43 


‘20.2 




4 




Mtl 


if 




0 

tr 


0 
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0 0 
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0 


im 


-iSg.4 . 




ft' 

1. J 




0 

. 



(II. C. 4-19) 



where the 



The state 



state 




variables are; 

X 2 ’ pitch angular rate state variables 
(section II. C. 4) 

^ 2 - Achieved Body Acceleration state 
variables (section II. C. 4) 

: output of Acceleration filter 

: input command in actuator network 

: tail incidence angle 

y^, 6p of control law are shown in 



figure 2.6. 
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Executing the program of Appendix D, using an 
input step function which represents the "IGee Command” and 
trim angle of attack a^ =0, the time response plots of 
pitch normal acceleration n^f pitch angular rate q and 
pitch tail incidence 6 for circular airframe are obtained. 

£r 

Figure 2.7 shows the pitch acceleration response 
of the circular airframe due to a one Gee acceleration 
command, when the Missile Aerodynamic Models are linearized 
about zero angle of attack (a). The response has a 0.5 
seconds time constant, steady state error 0.005 and 2.5% 
overshoot, which are according to the requirements referred 
to in appendix C. 

Figures 2.8 and 2.9 show the required body 
angular rate and control surface deflection, to achieve the 
acceleration response. These responses matched those 
presented in [Ref 2]. 

5. Design Approach and Analysis For Elliptical 

Airframe 

The design technique used for the elliptical 
airframe, is the same as that of the circular airframe. 

The analysis of the Pitch channel for the elliptical 
airframe was based again, on the time response plots. 

For the above purposes, a CSMP program was 
written, as in the case of circular airframe, referred to 



in Appendix D. 
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Fig. 2.7 Pitch Normal Acceleration (n^) vs Time (t) 
Uncoupled Pitch Channel 

Circular airframe (IGee Command, a^ =0). 



45 



Q(RRO/SEC) 




Fig 2.8 Pitch Angular Rate (q) vs Time (t) 

Uncoupled Pitch Channel 

Circular airframe (IGee Command, a =0) 

e 
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Fig. 2,9 Pitch Tail Incidence i Time (t) 



Uncoupled Pitch Channel 

Circular airframe (IGee Command, a^ =0) . 
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The equations of the uncoupled linear Pitch 
channel for elliptical airframe, which are used in the 
above CSMP program, are; 

a. Equation of Control Law 

The control law of the uncoupled linear Pitch 
channel for the elliptical airframe is very similar to the 
control law for the circular airframe, with only 
differences in the gain of acceleration and rate 
compensators. Elliptical control law is shown in figure 
2.6., taken from [Ref. 2]. 

(1) Acceleration Filter Equation . The 
Acceleration filter equation is the same as in section 
II. C. 4 



(2 ) Acceleration Compensator Equation 

Y= {-O.OOg/|-+lV (l.^053»l2,-x) (II. c. 5-1) 
utilizing, inverse Laplace Transformation 
and rearranging II. C. 5-1 it becomes: 

Vr -6 V- SSO.S^IMXIO^TIU (ii.c. 5-2) 

(3) Compensator Equation 

Spc _ -5-OT- 

v-x, ■ ^ 

(II.C. 5-3) 

Using inverse Laplace Transformation, it 



yields ; 

5pc - -3g3.75X^O^Cy-X^') -3.0y(Y-X<) (ii.c.5-4) 
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b. Actuator Equation (Figure 2.5) . 



As in section II.C.4,b 

c. Equation of Uncoupled Linear Pitch 
Aerodynamic Model 

The equation for pitch angular rate was 
derived, as in section II. C. 3, the only differences are in 
the^ constants W=2475 (lbs), 

lyy=790 (Slug-ft^) , C =0.015, C =-0.055, =0.02 and 

A ft > UK A Ji Jt A 

a 6 0 

P P 

=0.18, as referred in the Appendix A and Table I. 
a 

Manipulating according to the inverse Laplace 
Transformation rules and rearranging, one abtains. 

(3.66 S+0-6S8)/(-88-6MKi6^S^+ ^5.95X^6^S -+i) (II.C.5-5) 

Using state-space representation of a system, 
in which forcing function involves derivative terms [Ref. 

5; pp. 675-678.] It yields; 

= X 2 -Sp (II. c. 5-7) 

X 2 ~ 11.2Si5-X^+‘lT9.0^X'(o’x2-l5.'f91‘( <5’p (II.C.5-8) 
S = (II.C.5-8a) 

Also, Manipulating with same way the equation: 

, -2l2.>4ilXl5^S^->-gi.6 (II. C. 5-9) 

Jf> * -88.6qX<O^S^+15.e5X<O^S + 1 

Yielding 

2^ = h 528.3 Xld^iS'p 
22 = 11.2815-21 + 179.89 Xld^Z2 - ^87-35 -Sf> 

y\z - - ( ^1 + ^- ) 
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(II. C. 5-10) 
(II. C. 5-11) 
(II. C. 5-12) 



utilizing state-space representation, all of 
the above equation can be modeled in the following 



eighth-order system; 
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( 1 1. C. 5-13) 



where the 



The state 



state 



X 



1 ' 





variables are; 

X 2 ; pitch angular rate state variables 
(section II. C. 4) 

Y 2 ; Achieved Body Acceleration state 
variables (section II. C. 4) 

; output of Acceleration filter 

; input command in actuator network 
; tail incidence angle 

<5p of control law are shown in 



c 



Figure 2.6 
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utilizing, the above state variables system, a 
CSMP program was written (Appendix D) to achieve the time 
response plots of q and Sp. As input a step function 

is used, which represents a one gee command. The trim 
angle of attack (a^) is equal to zero. 

Figure 2.10 shows the pitch acceleration 
response of the elliptical airframe due to the command. The 
response has a time constant of about 0.5, steady state 
equal to zero, and overshoot o% , keeping the requirements 
referred to in Appendix C. 

Figures 2. Hand 2.12 show the required body 
angular rate and control surface deflection, to achieve the 
acceleration response. 

6. Comparison of Airframe Response 

Figures 2.7 and 2.10 show that the Pitch 
acceleration response of the circular and elliptical 
airframe are approximately the same having a 0.5 second 
time constant and negligible overshoots. 

Figures 2.8,2.9,2.11 and 2.12 show that to achieve 
the acceleration response the elliptical airframe requires 
less body angular rate and control surface deflection. The 
elliptical airframe requires less control surface 
incidences, because the transfer functions have a higher DC 
gain and poles which are located at a lower frequency due 
mainly to a more neutrally stable airframe. 
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Fig. 2.10 Pitch Normal Acceleration (n^) vs Time (t) 
Uncoupled Pitch Channel 
Elliptical airframe (IGee Command, 
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Fig, 2.11 Pitch angular rate (q) vs Time (t) 
Uncoupled Pitch Channel 

Elliptical airframe (IGee Command, a =0) 
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Fig. 2.12 Pitch tail incidence (5 ) vs Time (t) 
Uncoupled Pitch Channel" 

Elliptical airframe (IGee Command, a =0) 



54 



) ] is directly 



The DC t (^S/w) . ( Cjj * ^ ^in ~ 

a “6 a 5 

P P 

proportional to pitch control moment C and inversely 

P 

proportional to magnitude of stability margin in Pitch 

(i.e., C /Cjj or distance from center of pressure to 
™a a 

center of gravity) . Thus the nearly neutrally stable 

elliptical airframe is expected to have a higher gain than 

the stable circular airframe even though its control moment 

coefficient is slightly smaller. Higher DC gain will 

require less control surface incidence. The zeros of n^/'^p 

are directly proportional to C^ and the location of the 

ratio C / Cj, or the distance from the point of 
6 6 
P P 

action of tail forces to the center of gravity. 

Hence, for the elliptical airframe which has 

larger C^^ the zeros of are located at a higher 

frequency. Lower control surface incidences and higher 

frequency n /6 zeros of the elliptical airframe will 
^ P 

simplify the design of a rapidly responding Pitch 
autopilot. 

D. UNCOUPLED LINEAR YAW CHANNEL 

The purpose of the Yaw channel autopilot of a CBTT 
autopilot is to minimize sideslip angle or provide 
coordinated motion between Roll and Yaw channels. This was 
accomplished in two ways. First, the uncoupled Yaw channel 
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autopilot, (i.e. Roll and Pitch dynamic effect neglected) 
was designed as a regulator (i.e no guidance command and 
with rate and acceleration feedback) to help minimize 
sideslip angle. Second, to aid in sideslip control, the 
regulator was commanded from the Roll channel as explained 
in chapter III. 

1. Aerodynamic Model Control Law 

A block diagram of the uncoupled Yaw channel is 
shown in Figure 2.13, taken from [Ref. 2]. In this diagram 
the aerodynamic model and the Yaw control law are involved. 
Also, the Yaw control law without the accelerator is shown 
in Figure 2.14. These figures were abstracted from 
[Ref. 2]. 



The normal acceleration n is not used to 

^c 

command the CBTT autopilot. Instead, it is used for the 
design and analysis of the uncoupled channel. The command 
used by the coupled system is shown in dashed lines and is 
a yaw angular rate command, r 

The yaw control law is governed by missile body 
yaw angular rate ( r) and yaw normal acceleration n^. The 
yaw control law determines the required command 6^ 
to an actuator, which is approximated as a first order lag 
at 30 Hz (188.4 rad/ sec) . 

The aerodynamics, linearized about a trim a^ 
angle-of-attack equal to zero. 
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Uncoupled Yaw Channel 



CN 


IT) 


00 


CXD 


q 






CD 






LO 




CN 


CN 




d 


d 




CO 


IT) 






CO 




o> 


o> 




T— 


CO 




CO 


00 




d 


d 


0> 

E 


0) 


15 

o 


03 






k_ 


o 

o 


Q. 


u. 




— 


< 


d 


txi 



Q 

k_ 

O 

> 




CO 

'cn 



u 

> 



CN 

cn 

••H 

Cm 



58 



Yaw Control Laws 



2. Transfer Functions of Aerodynamic Model 

The aerodynamic characteristic of the Yaw channel 
are modeled by a second-order system, by the following 
equation, where the states are r, yaw angle and 6 , 
sideslip angle: 

Yr (q .^ol (lii) ' ( ^ + 57/3 C-n5y • 5 v ) (yocol /sec) 

^ - 575 

and also, by the equation: 

V\Y-(q^/W)« (Cy^-^ + 5T.3Cyjy.5y 1 (3«s) 

substituting equations II. D. 2-3 into 

last equation becomes; 

e - 5T.5 [-r + 2« (M (C^,.g + ST.3 Cy,,.iy)') J 

Taking the Laplace Transform of the equations, one 

obtains 

Ih 

gS-S1-3 + (c>,-e+57-3 (ySySv))) 

(II. D. 2-6) 

Rearranging above equations, in terms of q and o. , 
they become : 

5T.3-V+ (s- ^ 



C-ng • ^ 4 57.3 ♦ Jv ) 



(II. D. 2-5) 



(II. D. 2-1) 
(II. D. 2-2) 

(II. D. 2-3) 
II. D. 2-2, the 

(II. D. 2-4) 
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(II. D. 2-8) 




•Km 




Applying the Cramer's rule in the system of the 



equations II. D. 2-7 and II. D. 2-8 the transfer functions r/ 6 



and 6/6y are obtained. 






and 

±_ 

sp 



5’-( 






w In 

5=!-3 q <i 
l2i 



)• 



^ne 






52g3.2S-9^«^-^^(?v 

Ivv 






gr.5 q Cinoc 

lyv 



(II. D. 2-9) 



(II. D. 2-10) 



K 

q 

S 

d 

Izz = 



w 

C 



n , 



(57.3) . (32.3) 

V 

1650 (Ibs/ft^) 
3.14159 (ft^ 

2 (ft) 

810 (slug-ft^) 

2525 (lbs) 
-0.025 



= 0.48 



n 



6 = 

y 



8 = 



-0.050 

-0.065 



y« = 

y 



= 0 



0.021 

o 



As referred to in Appendix A and Table I. 
Substituting the aforementioned values into equation 
(II. D. 2-9), it becomes: 
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-(<^ec) (II. D. 2-11) 



r - 2.00S-5 

^Py ' -5il.54Xl6’5^- 2.8X10'^ 5 + i 
Introducing equation II. D. 2-10 into II. D. 2-3 and 

rearranging the last one, the transfer function of yaw 

normal acceleration in terms of the yaw control incidence, 

is obtained: 



Si 



( 5T-3Krqi?)^^ ^ ^ 3283.39(q^)^d ^ ^ 

_ [ ^ • fv# 



q j Gnct 



57 3 

w 






(II. D. 2-12) 

After substituting the values into the equation 
(II. D. 2-12), one obtains: 



S + IT.TC 



Jy ” -54.54X10^5 -ZSX^O-S + 1 



(gas/raot) 



(II. D. 2-13) 

The above transfer functions are used for finding 
the time responses of the yaw angular rate r and the yaw 
normal acceleration. 

Using inverse Laplace transformation, the equation 
II. D. 2-11 and II. D. 2-13 become: 
r +51.338X10^ f - J8.3335-r = -36.761- 5y - 2.75- 5y ( 1 1 . D. 2-1 4 ) 

My +5J.538X40^TiY-ia3335'rty= 2.387<5 y -325.629- (II. D. 2-15) 



The equations II. D. 2-14 and II. D. 2-15 involve 
derivative terras in the forcing function. Using rules of 
state-space representation of a system [Ref. 5: pp 
675-678] , they turn: 
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(II. D. 2-16) 



— X 2 - 36*76J*<5v 

^2 - 18.335-X^ - 51.338X‘10-3(2 ~S63.0-m X^O^•^y 



(II. D. 2-17) 



y - Xi 

it 3 ?2- 



(II. D. 2-18) 



(II. D. 2-19) 



= ^S.335-2l -51.33SX^0^i2- 2?J.S5?-Sy (II. D. 2-20) 



Odyr 2^ + 3.3?^(^y 



(II. D. 2-21) 



3. Equations of Yaw Control Law 

As mentioned before, the yaw control law is 

governed by missile body yaw angular rate (r) and yaw 

normal acceleration and it is used to determine the 

required command to the actuator. The yaw control law 

c 

at the flight condition of interest (i.e, Mach 3.95, 60kft) 
is shown in Figure 2.14, taken from [Ref. 2]. The rate 
compensation determines the high frequency attenuation and 
is used to minimize aerodynamic variations on the quality 
of regulation. The acceleration compensation determines 
the acceleration bandwidth via the time constant of the 
acceleration response due to a step command of acceleration 
at n^ . 



c 



a. Acceleration Compensator Equation; 




(II. D. 3-1) 



Inverse transforming and rearranging the 



II. D. 3-1, it turns into; 



V - -5Y +X.59^3 (ii.D. 3 - 2 ) 
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b. Rate Compensator Equation 

M8S (-fe-H) 






(Ytr) 



(II. D. 3-3) 



Utilizing inverse Laplace Transformation 
and introducing equation II. D. 2-16, the equation II. D. 3-3 
turns into: 

r 0.M85 (y + xO +4. 85 (Y + Xi) (II.D.3-3a) 

c. Actuator Equation: 

^ Syc (II.0.3-4) 

+ i 






Inverse transforming and rearranging the 
equation II. D. 3-4, it turns into: 

(Jy r - 188.^ • 5 y + JYc (ll.D.3-5) 

Using state-space representation, all the 
aerodynamic model and control laws can be modeled in the 
following seventh order system: 
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IIVc 



(II. D. 3-6) 

where the state variables are: 

X 2 : pitch angular rate state variables 

(section II. D. 2) 
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i' 




1 




The state X^/ 
Figure 2.14. 




yaw acceleration state variables 
(section II. D. 2) 

output of acceleration compensator 

input command in the actuator 
yaw tail incidence 
control law are shown in 



4 . Design Approach and Analysis of Circular Airframe 
The design technique used for the Yaw channel is 
the same as the technique used in the Pitch channel. The 
analysis of it, for circular airframe is based, on the time 
responses of the yaw angular rate r, the yaw normal 
acceleration and yaw tail incidence 6^, 

For purposes of analysis, a CSMP program was 
written (Appendix E) using the state-variable system which 
is represented by the block diagram of the yaw aerodynamic 
model and control law. 

Executing the referred program, using as input a 
step function, which represents the "one Gee Command" and 
trim angle of attack 3^=0/ the time response plots are 
obtained . 

Figure 2.15 shows the yaw acceleration response of 
the circular airframe. The response has a 0.39 seconds 
time constant, 7% overshoot and steady state error 0.018 
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Fig. 2.15 Yaw normal acceleration (n ) vs time (t) 
Uncoupled Yaw Channel ^ 

Elliptical airframe (IGee Command, a =0) 
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which are according to the requirements referred in 
Appendix C. 

Figures 2.16 and 2.17 show the required body yaw 
angular rate and control surface deflection, to achieve the 
acceleration response. Again these results matched those 
of [Ref. 2] . 

5. Aerodynamic Transfer Function for Elliptical 

Airframe 

The transfer functions for yaw angular rate and 
yaw normal acceleration for the elliptical airframe, were 
derived in the same way as in section II. D. 2. The 
only differences are in the following constants 



w 




2475 (lbs) 


Izz 


= 


853 (slug-ft 




= 


0.024 


Cn 


= 


-0.042 








"""a 


= 


-0.043 




= 


0.016 



As referred to in Appendix A and table I. 
a. Transfer function of yaw angular rate r. 

r . > 0.059, (ii.d.5-1) 

oV 59.858X10 .5^2.593X10-5 + 1 
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R(RflD/SEC) 

-0.039.000 0.035 0.050 0.075 0.100 0.135 O.ISO 0.175 0.300 0.335 0.2S0 0.375 0.300 







Tisec) 



Fig. 2.16 Yaw angular rate (r) vs Time (t) 
Uncoupled Yaw Channel 

Circular airframe (IGee Command, 3g“0) 
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Fig. 2.17 Yaw tail incidence (6^) vs Time (t) 
Uncoupled Yaw Channel 

Circular airframe (IGee Command, a^=0 ) 
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b. Transfer function of yaw normal acceleration n^ 



My _ »H.6BXH0-5^- Ti / j£u \ 

iy ' 59.83SX10'^-5^+2.593X<0^S+1 V rod ) 



(II. D. 5-2) 



Using inverse Laplace transformation, the equation 
II. D. 5-1 and II .D. 5-2 "becomes ; 

r+H3333X<6^f4^6.?ia'r = -29.213* Jy-9‘3<3.36X-|5-Jy (II -D. 5-3) 

-^43.333Xl6^7iy 4<6.?Jl? 'Wy = 1.<^15?-<Jy-118.653-^y (il.D.5-4) 

utilizing state representation of a system in 
which the forcing functions involve derivative terms, as in 
section II. D. 2, the above equations yield; 

Xi - X2- 29.221 ‘ J'y 
5(2 - -ium -(43.333X10^X2 
rr Xl 

32- '?3.015X>10^(5‘y 

?2 - “lG).^113*2i -43.333X10?22-'150*6633*Sy 
2^-H.951^*^y 
The aerodynamic model for elliptical airframe is 
shown in Figure 2.13. 

6. Yaw Control Equation for Elliptical Airframe 

The control law of uncoupled linear Yaw channel for 
elliptical airframe, is similar to the control law for 
circular airframe. There are differences in the gains of 
the acceleration and rate compensators, as is shown in 
Figure 2.14. 



{II. D. 5-5) 
(II.D.5-6) 
(II. D. 5-7) 
(II. D. 5-8) 
(II. D. 5-9) 
(II. D. 5-10) 
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(II. D. 6-1) 



a. Acceleration Compensator Equation 

V= ( 0 . 8393 /o.a 5 S+l) • 

Utilizing inverse laplace transformation and 
rearranging the equation II. D. 6-1, it turns into: 

V = -M V + i35?M (II. D. 6-2) 

b. Rate Compensator Equation 

(y-y) 

^ (II. D. 6-3) 

Inverse Laplace transforming and introducing 
the equation (II. D. 5-7) into the equation (II. D. 6-3), it 
turns into: 

« 

S>/^ = 0.608 (y+X<) +6.08(Y+XA (II. d. 6-4) 

c. Actuator Equation 

Same as the equation II. D. 3-5. 



Using state-space representation all equations of the 
aerodynamic model and control law, can be modeled in the 
follow seventh-order system: 
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(II. D. 6-5) 
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where the state variables are; 



Xj, Xj: 




The state » Y, 
2.14. 



yaw angular rate state variables 
(section II. D. 5) 
yaw normal Acceleration state 
variables (section II. D. 5) 

output of acceleration compensator network 
input command in the actuator 
yaw tail incidence 

•Sy of control law are shown in Figure 
c 



7. Design Approach and Analysis for Elliptical 

Airframe 

The design technique and analysis for elliptical 
airframe is the same as that, in section II. D. 4. Again, 
for purposes of analysis, a CSMP program was written 
(Appendix E) using equations of aerodynamic model and 
control law. 

Figure 2.18 shows the yaw acceleration response of 
the elliptical airframe. The response has 0.35 second time 
constant, overshoot about 6% and steady state error 0.01 
which are according to the requirerments , (Appendix C) . 

Figures 2.19 and 2.20 show the required body yaw 
angular rate and control surface deflection, to achieve the 
acceleration response. 
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NY (GEES) 




Fig. 2.18 Yaw Acceleration (n^) vs Time (t) 
Uncoupled Yaw channel 

Elliptical airframe (IGee Command, a^ 
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Fig. 2.19 Yaw Angular Rate (r) vs Time (t) 
Uncoupled Yaw Channel 

Elliptical airframe (IGee Command, a^=0) 



73 



DY(RRD) 




Fig. 2.20 Yaw Tail Incidence (6 ) vs Time (t) 
Uncoupled Yaw Channel^ 

Elliptical airframe (IGee Command, ^^=0) 
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1 



1 



E. 



LINEAR UNCOUPLED ROLL CHANNEL 



Generally, the purpose of a roll channel autopilot is 
to command the t signal, in order to roll the missile 
through an angle <t>, measured from the vertical axis. 

The roll channel of a coordinated bank-to-turn 
autopilot is commanded to roll the missile so as to put the 
preferred maneuver direction of the missile in the 
direction of the guidance acceleration command, while the 
pitch channel acceleration is commanded to produce the 
total magnitude of the guidance acceleration command. 

The desired maneuver plane acceleration should be 
attained as rapidly as the achieved body-fixed pitch 
acceleration. To accomplish this, the uncoupled roll 
channel autopilot was designed to have the roll angle time 
constant equal to the time constant of the normal 
acceleration achieved by the uncoupled pitch channel 
autopilot. This is not necessarily the optimal 
relationship betweem the two time constants. 

1. Roll Aerodynamic Model and Control Law 

A block diagram of the uncoupled Roll channel is 
shown in Figure 2. 21, taken from [Ref. 2]. In this diagram 
the aerodynamic model is involved. Also, the roll control 
law is shown in Figure 2.22. These figures were taken from 
[Ref. 2]. 

The roll control law is commanded by roll angle t 

is governed by roll angular rate p and roll angle <t> and 
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Roll Control Laws 



determines the required command 6^^ to an actuator which is 

c 

approximated as a first order lag at 30Hz (188.4 rad/sec) . 
The aerodynamics, linearized about a trim angle-of-attack 
(3g) equal to zero. 

2. Equations of Aerodynamic Model for Circular 
Airframe 

The aerodynamic model is represented by the 
following equations: 

pr (q^d/lxx)- iradkic^) (h-e.2-1) 

P(mdlsec) 

where : 

q = 1650 (Ibs/ft^) 
s = 3.14159 (ft^) 
d = 2 (ft) 

Ixx = 40 (slug-ft^) 

= 0.031 

As referred to in Appendix A and Table I. 
Substituting the above values in equation II. E. 2-1, it 
yields : 

P - S.03^G-(5'|2 (II. E. 2-3) 

3. Equations of Aerodynamic Model for Elliptical 
Airframe 

The aerodynamic model equations are similar to 
those for the circular airframe, with the only differences 
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0.023. 



in the constant I =110(slug-ft^) and C = 

X A « 

Hence they are: 

p - Cvadlsec^^ 

4) - P C'fad IseO 

4. Equations of Control Law for Circular 
Elliptical Airframe 

The yaw control law for circular and 
airframe is shown in Figure 2.22. 

a. Roll Angle Compensator Equation 



(II. E. 3-1) 
(II. E. 3-2) 
and 



elliptical 






‘ 2.2 



9 



+ i 






(II. E. 4-1) 



inverse Laplace transforming the equation II. E. 4-1 
it becomes; 

X = (II. E. 4-2) 

b. Rate Compensator Equation 






0.15065* 






(II. E. 4-3) 



utilizing inverse Laplace transformation, and 
minor manipulating the equation II. E. 4-3, one obtains: 

V r-5V+50.333XHO^(X-P) + 0.15S5(x-p) (ii.e.4-4) 



c. Psuedo Differentiator Equation 

0 . 0^303 

p 

4-+i 






(II. E. 4-5) 



Using inverse Laplace transformation and 
rearranging the equation II. E. 4-5, it becomes; 

Xf - “6X-) (II. E. 4-6) 
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d. 



Equation of Actuator Compensator 

This network of roll control law is different 



for circular and elliptical airframe. 

(1) Circular Airframe 

(II-E.4-7, 

Using, inverse Laplace transformation and 
rearranging the equation II. E. 4-7, it yields: 

4c ■ + (Y-X0 + ^5(V-X') (ii.e.4-8) 

(2) Elliptical Airframe 

Inverse Laplace transforming and 
rearranging the above equation, it becomes; 



= -isipc +o.5«993 (y->:))+6s.536(y-X)> 



(II. E. 4-10) 



e. Actuator Equation 

= (II. E. 4-11) 

4 i 

188.4 ^ 

Utilizing inverse Laplace transformation, the 
equation II. E. 4-11 becomes; 

•^^ 95-32 ( j j g 4_2 2) 

utilizing state-space representation, all the 
aforementioned equation of the aerodynamic models and 
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control laws for the circular and elliptical airframes, can 
be modelled in the state variables systems given below: 

For the circular airframe 
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(2) For the elliptical airframe 
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where the state variables are: 

P : roll angular rate 
<i> : roll angle 

X : output of roll angle compensator 
Y : output of rate compensator network 
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: output of pseudo differentiator 

: input command in the actuator 
c 

; roll tail incidence 

The state 6^ are shown in Figure 

c 

2 . 22 . 

5. Design Approach and Analysis of Uncoupled Linear 
Roll Channel for Circular and Elliptical Airframe 
The design technique and analysis is the same as 
that which was used for the Pitch and Yaw channel. 

For purposes of analysis, two CSMP program were 
written (Appendix F) using the state variable system 
represented by the aerodynamic model and control law. 

Block. diagram. 

Figures 2.23 shows the roll angle response of both 
elliptical and circular airframe. The time constant of the 
roll angle response is 0.55 seconds, the overshoot 3% and 
the steady state roll angle error equal to zero. 

Figure 2.24 shows the required body roll angular 
rate for both elliptical and circular airframes. Only the 
roll tail angular deflection figures 2.25 and 2.26 is 
different for the airframes, due to the method of 
compensating for a reduction in the elliptical aerodynamic 
roll gain. This is desirable for the aerodynamic roll gain 
(57.3qsdC^ /^xx^ large as possible to minimize 
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Fig. 2.23 Roll Angle ( 4») vs Time (t) 

Uncoupled Roll Channel 

Circular or Elliptical airframe 

(1 RAD Command, a =0) 

e 
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P(RflD/SEC) 







TISEC) 

Fig. 2.24 Roll Angular Rate (p) vs Time (t) 
Uncoupled Roll Channel 
Circular or Elliptical airframe 
(1 RAD Command. 3^=0) 
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T(SEC) 

Fig. 2.25 Roll tail incidence vs Time (t) 

Uncoupled Roll Channel 

Circular airframe (1 RAD Command, 

a^=0) 
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Fig. 2.26 Roll tail incidence vs Time (t) 

Uncoupled Roll Channel 
Elliptical airframe (1 RAD Command, 
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control surface motion. It is also desirable to have as 



large C as possible to minimize the effects of 
aerodynamic control cross coupling. The circular airframe 
has a considerably larger aerodynamic roll gain due to a 
much smaller roll inertia and a larger control derivative 
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III. LINEAR DESIGN AND ANALYSIS OF COUPLED AUTOPILOTS FOR 



CIRCULAR AND ELLIPTICAL AIRFRAMES 



A. GENERAL 

The linear design and analysis technique began with 
uncoupled autopilot channels. The uncoupled autopilot 
design technique was classical, using a combination of 
frequency response and root locus techniques [Ref. 2], to 
achieve practical bandwidths (i.e, sufficient high 
frequency attenuation) and in turn provides the range of 
required missile body angular rates and control motions, as 
mentioned in section II. 

In this section, CBTT control laws are analyzed to add 
control coupling for coordinated missile motion. A measure 
of sideslip control is obtained by applying a roll angle 
command to the linearized CBTT autopilot. The relative 
stability of the autopilot branches and means for improving 
stability are discussed. An examination of the autopilot 
sensitivity to aerodynamic cross-coupling is made. 

B. AERODYNAMIC MODELS FOR THE CIRCULAR AIRFRAME 

The aerodynamic models of pitch and roll-yaw channels 
for the circular airframe are shown in Figures 3.1 and 3.2, 
taken from [Ref. 2]. The models are linearized and were 
developed for stability studies in the frequency domain. 



83 






' 4 ^ 




The assumptions for linearization are the same, as those in 
section II.C.l. 

Pitch and Lateral channels are coupled via constant 
missile role rate P^. Also, the linear model has inertial, 
kinematic and aerodynamic cross couplings, which will be 
discussed further, in the non-linear model. 

1. Pitch Channel 

The linear aerodynamic model of the pitch 
channel consisted of the uncoupled model mentioned in 
chapter II, which is coupled with the Roll-Yaw channels via 
constant roll rate (Pg) » as shown in F igure 3.1. 

In this work, the Pe was set to zero, because 
when Pe=0, the Qe (i.e equilibrium pitch rate) has been 
found to have negligible influence in the lateral model 
(i.e roll-yaw) . 

The equations, which represent the pitch 
aerodynamic model given in (Ref. 2] are the same as the 
equations of the uncoupled linear model, mentioned in 
chapter II. 

2. Lateral Channels 

The linear lateral aerodynamic model consisted 

of the uncoupled roll and yaw aerodynamic channel, coupled 

via the linearized aerodynamic derivatives C , C and 

S ^ 6 ^ 

. Also, it is coupled with the pitch dynamic model via 
^R 
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the Pe, as previously mentioned. The lateral dynamic model 
is shown in Figure 3.2, as given in [Ref. 2]. 

To find the equations which represent the 
lateral dynamic model, Pe and Qe are set to zero because of 
their negligible influence in the lateral channel. 

The aforementioned equations are as follows: 







(II.B.2-1) 




z 6 * Cy^ 4- Cyjy- <Fy (<ie^) 


(III. B. 2-2) 


where : 




T r C^ieglscc^) 


(III. B. 2-3) 




£2 - + C-Ti^y ' Ei 


(III. B. 2-4) 




P - 5^3 9 tq /Ixx (ae^/sec^) 


(III. B. 2-5) 




P/si.3 (rotdUeO 


(III.B. 2-6) 




E3r (e^y • Sy + Sfi Cdeg) 


(III. B. 2-7) 




Eqz £ 5 f Ce^.^ (deg) 


(III.B. 2-8) 




6"= BA/5f*3-T+K-(9;^Edv\/) (deglsec) 


(III.B. 2-9) 


where : 


k 


= 0.48 




S 


= if(ft)^ 




d 


= 2 . 0 (ft) 




q 


= 1650 (Ibf/ft^) 




w 


= 2525 (Ibf) 




Ixx 


= 40.0 (slug-ft^) 




Izz 


= 810.0 (slug-ft^) 




Cy 

V 


= 0. 22 
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-0.082 



= 0.018 
C 

^6 = -0.019 



= - 0.053 

= -0.016 

e = -0.009 

= 0.035 
A = 10° 



As referred to in Appendix A and Table I. 

Substituting the values into the above equations. 



they yield: 

T r -13.95M -<5'p- 3S.?T-JY(de^litc^J{iii.B.2-10) 

p I - 133.66 •<?+5^g.7?g-^|^ -23TC^8iy(<l«3l3€C^) 2-11) 

g :r0.0l^M5 p-r-0.080g 5 *l' 0 .:iifSY(deg/$cc) (ih-b.2-12) 

4>= P/ST.3 (tadlscO (III. B. 2-13) 



C. CBTT AUTOPILOT CONTROL LAWS FOR CIRCULAR AIRFRAME 
The control laws, which were used by the CBTT 
autopilots of the circular and elliptical airframe, are 
shown in Figures 2.6 and 3.3, as given in [Ref. 2]. 

The pitch control law for the circular airframe is the 
same as determined in the uncoupled pitch channel study 
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Actuator 
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CBTT lateral control law 



section II. B. 2, and therefore, the equations which 
represent it, are the same as in section II. B. 3., except 
the equation for pitch tail incidence angle which is: 

cfp z -18H.4-(Jp4'IO^Q5.3 (deg/sec) (iii.c.i) 

The uncoupled roll and yaw laws are coupled via a 
cross-coupling branch, shown in solid line in Figure 3.3, 
in order to consist the linear coupled lateral control law. 

In the roll actuator of the lateral control law, the 
frequency has been changed from 110 rad/sec to 60 rad/sec, 
for improving of autopilot stability. 

The autopilot cross-coupling branch between Roll and 
Yaw channels has been added to provide coordinated motion. 
Coordinated motion or zero sideslip angle is achieved by 
directing the body fixed pitch axis of the missile of the 
missile velocity vector, so that there is no component of 
missile velocity along the body fixed yaw axis of the 
missile. Figure 3.4, as given in [Ref. 2] shows the 
attitude of the missile body with respect to its velocity 
vector V. 

When commanding an upward maneuver (i.e, * =0), the 
missile body moves upward with its pitch axis directed at 
the velocity vector until it reaches the desired maneuver 
level or angle-of-attack. No roll motion is required to 
maintain coordination for this maneuver. 
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0=0 



0 = 45' 




Fig. 3.4 Coordinated Missile Motion for 
Coordinated bank-to-turn (CBTT) 
and limited bank-to-turn (LETT) 
Control policies 
96 



Foe maneuvers in the 4=45 or 90 degree directions, the 
missile body rolls the desired 4 and moves its pitch axis 
directed at the velocity vector, until it achieves the 
desired maneuver level. For maneuvers with 4 higher than 
90°, the missile must roll about the velocity vector, while 
the Yaw channel directs the pitch axis towards the velocity 
vector for minimum sideslip angle. 

The missile is forced to roll about its velocity 
vector, when the desired maneuver direction is in the 
negative angle-of-attack direction, because it is desired 
to avoid negative angles-of-attack. These maneuvers are 
illustrated in Figure 3.4. 

The coordinating command r^,, is a yaw angular rate from 
a rate gyro. The command is equal to - Pa^, where 3 g=a 
(Figure 3.3) in the coordinating branch is the equivalent 
angle-of-attack, which is exactly equal to angle of attack 
(a) , in the linear studies. The equations which represent 
the lateral channel are as follows: 

1. Roll Channel 



a. Roll Command Compensator Equation 



Same equation as in section II. E 



Xr (ct)c-4)) 



(III. C. 1-1) 



b. Roll Rate Compensator Equation 




. 4 



s 



(x- P/si s) 



(III. C. 1-2) 
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Inverse Laplace transforming and rearranging 
the equation (III. C. 1-2), it becomes; 

--51,+ 0.0 505 ()(-f>/5r-3) + 0.T55 (X- P/ST-S) 

c. Psuedo Differentiator Equation 

Xl- O-Oii S ( p/ 5 ^ 3 ^ 



(III.C. 1-3) 



(III. C. 1-4) 



+ i 



Utilizing inverse Laplace transformation and 
minor manipulating the above equation, one has: 

r - 0 + 0.07? ( P/513) (III.C. 1-5) 

d. Actuator's Compensator Equation 

The frequency of the compensator was changed 
from 110 rad/sec to 60 rad/sec for purposes of stability, 
as referred in section III.C. The equation represents the 
network is: 






60 



+ 1 



1 5 



4 1 






(III.C. 1-6) 



Using inverse Laplace transformation, the 
equation III.C. 1-6, one obtains; 

Sf. = -15 + 0.25 (V-X,) +15 (Y-X,) liii.c.i-7) 

e. Actuator Equation 

4c (III.C. 1-8) 

- +• i 

utilizing inverse Laplace transformation and 
rearranging the equation III.C. 1-8, it yields. 
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M "V AO+^5- 2)2 ^Rc ( 1 1 1 . c . 1—9 

All the referred to above state variables/ are 



shown in figure 3.3. 

2. Yaw Channel 

a. Acceleration Compensator Equation 



0.52 

0.2 5+ 1 



'My 



(III. C. 2-1) 



Inverse Laplace transforming, equation 
III. C. 2-1 becomes: 

V 3 -5Y +1.G lly (III. C. 2-2) 

b. Actuator Compensator Equation 

^/c- (y-0.M58a-P/5T.3+t) (III. C. 2-3) 

Utilizing inverse Laplace transformation and 
rearranging the equation III. C. 2-3, it yields: 

- O.ijSsCv+ylsT.S-O.MSS a. F>/5f3)+ M.3S(y+ r/5T.5-C.i(58.a.P/5f3') 

(III. C. 2-4) 



c. Actuator Equation 

<5^yr +^0?9S.3•5Yc ( III. C. 2-5) 

Making use the state-variable representation, 
the above equations, of the aerodynamic models and control 
laws, can be modelled in the following, state variable 
equation system: 
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(III. C. 2-6) 



where the state variables are: 



r 

P 

4> 

3 

X 



R 



yaw angular rate 
roll angular rate 
roll angle 
sideslip angle 

output of roll command compensator 
output of roll rate compensator 
output of pseudo differentiator 
input command in roll actuator 



: roll tail incidence 

Y : output of yaw acceleration compensator 

Sy ‘ input command in yaw actuator 
c 

6y ; yaw tail incidence angle 

All the referred before state variables are 
shown in Figures 3.2 and 3.3, 
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D. ANALYSIS OF LINEAR CBTT AUTOPILOTS FOR CIRCULAR 

AIRFRAME 

In the analysis of the linear CBTT autopilots it is 
assumed that the missile is initially in the desired 
maneuver plane, is trimmed at ten degrees angles-of-attack, 
the equilibrium roll rate Pe and equilibrium pitch rate Qe 
were set equal to zero. 

For purposes of analysis, a CSMP program was written 
(Appendix G) using the equations of the aerodynamic models 
and control laws of the linear CBTT autopilots for the 
circular airframe. 

Figure 3.5 shows the sideslip angle for the linearized 
CBTT autopilot, having a maximum value of 2.32°. Since the 
model is linear, and therefore the magnitude of the 
sideslip is directly proportional to the magnitude of the 
input command, a one radian roll angle command was used for 
convenience . 

The coupled autopilots have a roll angle response which 
is the same, as the uncoupled roll channel response, 
differing slightly in the overshoot, as shown in Figure 

3.6. 

The aforementioned values and plots are exactly the 
same as the plots given in the [Ref. 2] therefore its 
results are verified. 
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BETA (DEG) 

0.2S 0.00 0.2S O.SO 0.7S 1.00 I.2S I.SO I.7S 2.00 2.2S 2.S0 




Fig. 3.5 Sideslip angle ( 5 ) vs Time (t) 

Linear CBTT Autopilot (1 RAD Roll 
Command); Circular airframe (a =10 deg, 
Pe=Qe=0 deg) ® 
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PHURflD) 




Fig. 3.6 Roll angle {^) vs Time (t) 

Linear CBTT Autopilot (1 RAD Roll 
Command) Circular airframe (ag=10 
deg, Pe=Qg=0 deg) 
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Comparing the state variable form system with these 
previous discussed in the section II. D and II. E, this is 
the twelfth order system because it includes roll and yaw 
autopilots and the matrices A,B are different than the A,B 
matrices at the uncoupled roll and yaw autopilot, because 
of the existence of the coordination branch. 

E. AERODYNAMIC MODELS FOR ELLIPTICAL AIRFRAME 

The dynamic models for elliptical airframe are the same 
as those of circular airframe and are shown in Figures 3.1 
and 3.2. 

The differences of the circular and elliptical models 
are in the following constants and linearized aerodynamics 
derivatives : 

2475 (Ibf) 

110 (slug-ft^) 

790 (slug-ft^) 

0.015 

-0.054 

0. 014 

0.024 

-0.039 



w = 

Ixx = 
Izz = 
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- 0.010 



3 = -0.027 

Cl 

6 ^ = 0.023 



as referred to in Appendix A and Table I. The equation 



which represent the aerodynamic model, are as follows: 



r Z i?.OM6<3- 6 29.3262 -<5v 
p r -•i95.'?^0^-5+^2M.209f5g-5M.OH*5y 
^ - P)57-3 

^ r 0.n452 p- V + 0.05428T •^*O.Of5C7S’'(Jy 



(III.E. 1) 
(III.E.2) 
(III.E. 3) 
(III.E. 4) 



F. CBTT AUTOPILOTS CONTROL LAWS FOR ELLIPTICAL AIRFRAME 
The control laws, which were used by the CBTT 
autopilots of the elliptical airframe are the same as those 
of the circular airframe and are shown in Figure 2.6 and 

3.3. 

The differences which exist in the control laws of the 
two airframes are the gain Kyp of the coordinated branch, 
and the gains of the actuators. 

Hence, the equations that represent the control laws 

are ; 

1 . Pitch Channel 

The linear pitch control law equations for 
the elliptical airframe, are the equations of the Pitch 
uncoupled channel, as given in section II. C. 5. 
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2. Roll Channel 



Same equations as in par, III.C. 

a. -X =-8X+ 17.6 (4>c-4>') (II. C. 1-1) 

b. Y, r -5 0.0505 (X-P|s?.S)-)-0.?55 (x-P/57.3) (II. C. 1-3) 

c. + O.OTS (P/ 57 . 3 ) 

d. Actuator Compensator 

T? ^ 

Utilizing inverse Laplace transformation and 

rearranging the equation III. F. 2-1, it becomes: 

« 

= -i5^R^,+ 1.0^25 (VpXf) 4 62.53 (Y<-V<) (iii.f. 2-2) 

e. i'R r -.182.4 + 10. 7053 (III.C. 3-8) 

3. Yaw Channel 

a. Acceleration Compensator Equation 



(II. C. 3-5) 
(III.F. 2-1) 



y= 



0.839 






(III.F. 3-1) 



+ 1 



Inverse Laplace transforming and rearranging 
the equation III.F, 3-1, it yields: 

Y: -9Y +5.356 >ly (III.F. 3-2) 



b. Actuator Compensator Equation 

^ _ai08£+_6j^ ( Y+ W5T.3 -dP/Sfs) 

5 

(III.F. 3-3) 



Utiliing inverse Laplace transformation and 
rearranging the above equation, one obtains: 

<5y^ : 0.608 (Y4^-|L) +6.08 (Y+^3-^j ) (III.F.3-4; 
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c. Actuator Equation 

Same as in par. III.F.2.e 

Utilizing the state-variable representation, 
the equation which represent the aerodynamic model and the 
control laws, can be modelled in the following 
twelfth-order system; 
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(III.F. 3-5) 



where the state variables are: 
r : yaw angular rate 

p : roll angular rate 

"t* : roll angle 

S ; side slip angle 

X : output of roll command compensator 

: output of roll rate compensator 

: output or psuedo-dif f erentiator network 

: input command in roll actuator 



roll tail incidence angle 
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Y 



output of acceleration compensator in yaw 
control law 



6 



Y ' 

c 



6 




input command in yaw actuator 
yaw tail incidence angle 



All the state variables are shown in Figure 3.2 and 
3.3. The above system is different in the matrix A and B 
from the system discussed in chapter II for uncoupled Roll 
and Yaw autopilots, because of the coordination branch 
equation . 



G. ANALYSIS OF LINEAR CBTT AUTOPILOTS FOR ELLIPTICAL 

AIRFRAME 

In the analysis of the linear CBTT for elliptical 
airframe, the same assumptions, which were made for the 
circular airframe, are valid. 

A CSMP program was written (Appendix G) with the above 
equations, which represent the aerodynamic models and the 
control laws of the linear CBTT autopilots. 

Figure 3.7 shows the sideslip angle B for the 
linearized CBTT autopilot, having maximum values of -0.46^ 
and 0.26°. 

Figure 3.8 shows the roll angle response which has an 
overshoot of 4 percent and is about the same with the roll 
angle of uncoupled roll autopilot. The above plots verify 

the results, given in the [Ref. 2] . 
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Fig. 3.7 Sideslip angle (B) Ys Time (t) 

Linear CBTT Autopilot (1 RAD Roll 
command) Elliptical airframe (a =10°, 
Pe=Qg=0°) 
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Fig. 3.8 Roll angle ( 4 .) vs Time (t) 

Linear CBTT Autopilot (1 RAD Roll 
command) Elliptical airframe 

(ag=10°, Pe=Qg=0°) 
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H. RELATIVE STABILITY OF CBTT AUTOPILOTS OF CIRCULAR AND 



ELLIPTICAL AIRFRAMES 

When the coordinating branch gain Kyp was set to zero 
in CBTT autopilot of the circular airframe, the maximum 
sideslip angle increased to 2.8 degrees. Therefore, for 
the circular airframe the coordinating branch is not very 
effective in helping the yaw autopilot to reduce sideslip 
angle . 

If kyp is increased to unity, as it is for the CBTT 
autopilot of the elliptical airframe, the maximum sideslip 
remains below one degree for 3 seconds but the autopilot is 
unstable. Modifying of the circular airframe will be done 
in the non-linear CBTT, for minimization of the sideslip 
angle . 

Removal of the aerodynamic cross-coupling (i.e C^^ , 

C^ , Cj^ ) had no effect on the sideslip angle. This 

showed that the aerodynamic cross-coupling plays an 
indirect role in determining the quality of sideslip 
control by determining the relative stability of the 
coordination branch or the magnitude of the coordination 
gain K^p. The magnitude of sideslip angle is dependent on 
the nulling effects of two parallel paths shown in Figure 
3.9, as given in [Ref. 2]. 
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The contribution of yaw acceleration n^ to the maximum 
0 is negligible and therefore neglected. S is formed 
mainly by the subtraction of the kinematic paths of a^p and 
the yaw angular rate r. 

The coordination is obtained in the CBTT control law by 
commanding the yaw autopilot with a yaw angular rate 
command r^ (Figure 3.2) of Kyp a^p, which forces r to be 
equal to a^.p and therefore nulling as shown in Figure 
3.9. The nulling process will be accomplished more 
efficiently if Kyp =1.0. The reason the sideslip is not 
nulled completely is that n^ is not zero and r can not 
equal r^ over all frequencies. 

The linearized CBTT autopilot of the circular airframe 
was unstable, when the uncoupled channel control laws used 
with a coordinating branch gain Kyp of unity. The 
autopilot was stabilized by decreasing to 0.458. The 
actuator of roll control channel has a frequency of 60 
rad/ sec . 

The stability of the CBTT autopilot for the elliptical 
airframe is considerably better, due to the setting of Kyp 
at unity, which results in a substantial decrease in the 
magnitude of sideslip angle as shown in Figure 3.9. 
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I. CONCLUSIONS 

1. The linear coupled autopilots verify the results of 
the uncoupled autopilots and keep the sideslip angle small 
for coordination motion, (section II, D. and II. E) 

2. Further Conclusions about the relative stability of 
the autopilots, about the gain K^p of the coordination 
branch and about the sensitivity of the autopilots to the 
aerodynamic cross-coupling are given in chapter 7 of [Ref. 
2 ] . 
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IV. NONLINEAR ANALYSIS OF CBTT AUTOPILOT FOR CIRCULAR AND 



ELLIPTICAL AIRFRAME 



A. GENERAL 

The three dimensional nonlinear aerodynamic models 
used for the analysis are shown in the Figures 4.1 and 4.2. 
For the configuration under consideration, the assumptions, 
which have been found to be consistent are referred to in 
chapter 5 of [Ref. 2]. 

Nonlinear aerodynamics and mass parameters values are 
presented in Appendixes A and B, as taken from the [Ref. 

2]. The same flight conditions (i.e 60kft altitude and 
M=3.95) as in linear CBTT. 

In this section the control laws for the non-linear 
CBTT autopilots are discussed. The results for commanding 
the CBTT autopilots for both airframes are shown and the 
desired aerodynamic model to enhance CBTT performance is 
determined. The performance of the CBTT autopilot for the 
elliptical airframe is determined, when no lateral 
aerodynamic cross-coupling exists and with ideal airframe 
dynamics. Also the performance of the ?qCBTT autopilot for 
the circular airframe is determined, when no inertial and 
kinematic cross-coupling exist into the pitch autopilot. 

The conclusion of the nonlinear analysis are stated. 
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ig . 4.2 Nonlinear Lateral (Roll/Yaw) channel dynamic 
model 
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B. CONTROL LAWS FOR ELLIPTICAL AIRFRAME 

The control laws used for the following nonlinear 3-D 
studies were the same as those used for the linear studies 
in section III, which are shown in Figures 4.3 and 4.4 
taken from [Ref. 2 ] , except for the gain as shown in the 
bold line of the coordination branch in Figure 4.4. The 
new gain a is held constant at one degree magnitude for 
angles of attack less than one degree positive and greater 
than negative 5 degrees. For angles-of-attack greater than 
one degree positive, the gain a is equal to the angle of 
attack. This maintains coordination for very small 
angles-of-attack. 

Gravity effects were not included in the linear 
studies, because it was assumed to have a negligible 
influence on autopilot stability and response perturbations 
about a missile trim conditon. 

Gravity effects were included in the nonlinear 
studies, where the missile body-fixed yaw axis will be 
subjected to the full force of gravity and may therefore 
have a significant influence on sideslip. 

1 . Pitch Control Law 

The pitch control law for the elliptical airframe 
is the same that is used in the linear studies in section 
III.F., except in the command applied to the channel. 
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CBTT Nonlinear lateral control 




In inertial rectangular coordinates, the 

acceleration commands, are: 

a. n^ (acceleration command in inertial Zy 

direction as in figure 2.4) - ~ cos^ IV. B. 1-1 

c 

where; r \2 • guidance command (gees) 

c 

-cos<l> : anti-gravity bias command (gees) 

9 : Elevation Euler angle = ( qcos4'-rsin‘l>) dt 



b. ny acceleration command in inertial 
c 

direction, which is equal to zero, because it is desired 
coordinated motion. In the polar coordinates, the 
acceleration commands are: 



a. 'Vic Z - [ ( 



(IV. B. 1-2) 



b. 




(IV. B. 1-3) 



hence = nn and n =0,1,2... 

Since the pitch control law of the elliptical 
airframe does not have an integration in the acceleration 
error path (i.e in the acceleration compensator as in the 
case of circular airframe) and requires a gain in series 
with the acceleration command shown in Figures 4.3 it was 
necessary to modify the antigravity command as follows to 
assure an anti-gravity bias of just one gee. 
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n (elliptical airframe inertial acceleration command) - 
c 

= n^ - 0.913 cose (III.B.1.4) 

c 

The autopilot was tested, using a commanded 2 gee 
(0°, 180°), that denotes an inertial command of 2 gees, 
which is first applied in the 0° or upward direction at = 
2 seconds for 3 seconds and at the fifth second is applied 
in the 180° or downward direction. Since at first both the 
missile roll angle and roll angle command are at zero 
degrees, there is no roll motion and the missile turns 
upward as a skid-to-turn controlled missile. At fifth 
second the missile is commanded to roll through 180 degrees 
while moving in a coordinated manner in yaw and roll to 
minimize sideslip angle and prevent or minimize negative 
angle-of-attack . 

The other equations, which represent the nonlinear 
pitch control law are the same as for linear studies 
section III.C. and the state variables are shown in Figure 

4.3. 



The equations are: 

X z -150 X + 150 

<f r -6V- 0.5305 \ + O-^gX 

z -0.383T(Y+4l5T-3)-3.07fY-^l5?.3) 

Sf ' *5p + 



(IV. B. 1-1) 
(IV. B. 1-2) 
(IV. B. 1-3) 
(IV. B. 1-4) 



122 




m 




where : 



n„ : achieved body-fixed acceleration in Z axis 

n^ : inertial acceleration command, 

c 



2 . Lateral Control Laws 

The lateral control laws are the same as those 
used for the linear studies, except the modification in the 
coordination branch. 

The command ^ for the roll control law is zero 
c 

for the first 5 seconds and 180° affterwards, as derived in 



IV.B.l. 



The equations which represent the lateral control 



laws are the same those of section III.F. 

0 

r -L| V 2 -t- 3.556 >tYe (iv.b.2-1) 

Si,-- 0.608 + ry, 

^ ^ 5T-3 . (5T.3r ^ V 2 5T.3 C5T-3)* / 

z 0.608 fy2+x__5sp ivP 

ST.3 3383.23 32S3. 






(IV. B. 2-2) 



where : 

\ 

Ot = a J 

and 

a<i° 

at= 0 J 

1, = -8X(-1«<t>+JW't>c 

<jy = -138.4 <Ty +10W5.3 S'lc 

% Z -5lf(y-0.05O33(i:,-p)B.3)+air5s (x,-p/5t-3) 

XjZ-6Xi + 0.048 (p/5«') 

41.04i5 CyH“X^3+62.58 (yt-Xi) 



(IV. B. 2-3) 
(IV. B. 2-4) 
(IV. B. 2-5) 
(IV. B. 2-6) 
(IV. B. 2-7) 



r -igg.il-Jp 4 ^0'^95.3 -<5ec 




The state variables are shown in Figure 4.4. For 



purposes of analysis, a CSMP program was written (Appendix 
H) . 

C. AERODYNAMIC MODELS FOR ELLIPTICAL AIRFRAME 

The aerodynamic models of Pitch and Roll-Yaw channels 
for the elliptical airframe are shown in Figures 4. 1,4. 2. 
The airframe dynamics under consideration, have a level of 
complexity sufficient to determine the critical areas of 
concern regarding the stability and response of CBTT 
control . 

1. Pitch Channel 



channel consisted of the uncoupled model, coupled with the 
Roll-Yaw channels via the inertial cross-coupling of rp 
into q and the kinematic cross-coupling of -S.P into a. 

Also the antigravity command cos4 cos0 is inserted into n . 
This model is shown in Figure 4.1. 



The nonlinear aerodynamic model of the Pitch 



The above terms consist of the nonlinear terms of 



the pitch channel. 



The equations which represent the pitch 



aerodynamic model are: 



z -q Cn /w (g«es) 

q = Pr(5?.3+5^.3q;^d(H llyy 



(IV. C. 1-1) 



(IV. C. 1-2) 
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(IV. C. 1-3) 



cx + K.'vi^g- P.5/573 

where : 

1845 

K = = 0.48 

V 

q = 1650 (lbf/ft 2 ) 

S = IT (ft^) 

d =2 (ft) 

w = 2475 (Ibf) 

lyy = 790 (slug-ft^ 

Functions nonlinear 

functions shown in Appendix B that vary with a and 6p. 

Substituting the values into the above equations, 
they yield: 

•VJ2^= -J.3333 ♦ C/v( 

q r 0.17^52 * p-f f ^Sigs^J Cm 

<i - 9 -0.6399? Cm- 0.0^7^52 p. 5 

2 . Lateral Channels 

The nonlinear aerodynamic model consisted of the 
uncoupled roll and yaw channels, which coupled via 

B 

C and . Also the lateral channel coupled with 

the lateral channel coupled with the pitch channel via the 
inertial cross-coupling of -qp into r and the kinematic 
cross-coupling of ap into 8. The above term is the 
nonlinear part the lateral pitch channel. The nonlinear 
lateral dynamic model is shown in Figure 4.2. The 



(IV. C. 1-4) 
(IV. C. 1-5) 
(IV. C. 1-6) 
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equations that represent the aforementioned model are as 



follows : 

I ^ C$«es) (IV. C. 2-1) 

6 r K + a»p|sT.3 -f (cie^lsec) (iv.c. 2 - 2 ) 

r :: -q-PfST.^ +5T.5q^d C • ^ + 



P 

<i> 



5T.3 q ^ *^y + ) / Iax 

P/57.3 ('fouilsec) 

P (deglsec) 



(IV. c. 2-3) 
(IV. C. 2-4) 
(IV.C.2-4a) 
(IV. C. 2-5) 



where : 

K = 0.48 

q = 1650 (lbf/ft 2 ) 

S = u (ft^) 

d =2 (ft) 

w = 2475 (Ibf) 

^xx ” (slug-ft^ 

I = 853 (slug-ft^) 

zz 



as referred to in Appendix A. 



Functions C„ (a) , C„ (a) , C (a) , C (a) , C 
C. (a ) , C (a) and C (a) are nonlinear functions which 



vary with a angle-of-attack and referred in Appendix B. 

Substituting the values into the referred before 
equations one obtains: 



126 



1.5333 +Cyjy5y) 

S = 0. G3998 (%’6-¥ Cyjy .<^Y ) +0.01^^52 a-p- r 
f z -0.01?95 (^.P+696.U2 <5V +%p •S‘e) 

p r 5400.397? + Qj • Je) 



{IV. C. 2-5) 



(IV. C. 2-9) 



(IV. C. 2-7) 



(IV. C. 2-8) 



(IV. C. 2-6) 



D. ANALYSIS OF NONLINEAR CBTT AUTOPILOTS FOR ELLIPTICAL 
AIRFRAME 

For purposes of analysis, a CSMP program was written 
(Appendix H) with the aforementioned equations of the 
aerodynamics models and control laws of the nonlinear CBTT 
autopilots . 

To decrease computer time, appropriate initial 
conditions were inserted. The initial conditions are: 
a : angles of attack = 2.41 degrees 
6 p ; pitch tail incidence = 0.658 degrees 
9 : pitch Euler angle = 0.636 rad 

X : output of pitch acceleration feedback lag = Igee 

y : output of pitch acceleration compensator = -0.0105 

6 : pitch actuator command = 0.0148 rad 

^c 

The achieved maneuver plane acceleration n is 
calculated from the body-fixed accelerations n„ (Figure 
4.6) and (Figure 4.7) and the roll angle 41 (Figure 4.13) 
as follows: 
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Fig. 4.5 Achieved Inertial Acceleration (n^) vs Time ^t) 
CBTT of Elliptical airframe: 2 Gees (0 , 180 ) 
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Fig. 4.6 Achieved Body-fixed Acceleration (n„ ) vs 

B 

Time (t) CBTT of Elliptical airframe; 2 Gees 
(0°,180°) 
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Fig. 4.7 Achieved Body-fixed Acceleration (n„ ) vs 

B 

Time (t) CBTT of Elliptical airframe; 2 
Gees (0°,180°) 
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(IV.D-1) 



Figure 4.5 through 4.16 show the responses of n , 

Hy , a/ q# 9t tf '5 , 6 and 6j,. Also these figures 
and table II show the time constants and the overshoots of 
the referred before variables. 

Figure 4.5 shows that the achieved inertial 
acceleration n„ doesn't satisfy the requirement of the 
overshoot (i.e below 10 percent) and it is not acceptable. 
The results of this analysis verify the result of the [Ref. 
2 ]. 



E. CBTT PERFORMANCE WITH IDEAL AIRFRAME DYNAMICS FOR 

ELLIPTICAL AIRFRAME 

The purpose of the following simplifications to the 
airframe dynamics model is to isolate the critical 
cross-coupling paths, which have caused, the transients in 
the maneuver plane acceleration responses. 

The dynamic model without the coupling paths will be 
referred to as ideal dynamics. Although ideal dynamics are 
not physically attainable, it is a useful goal for both 
autopilot and airframe designers. 

The same guidance commands are applied to the CBTT 
autopilot of the elliptical airframe as in par. IV.B.l 
(i.e, 2 gees (0°, 180°) but with lateral aerodynamic 
cross-coupling removed. That means, the aerodynamic 
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Fig. 4.8 



Angle of Attack (a) vs Time ( t) CBTT of 
Elliptical airframe; 2 Gees (0°, 180°) 
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Fig. 4.9 Sideslip angle (g) vs Time (t) 

CBTT of Elliptical airframes; 2 Gees 
(0°, 180°) 
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Fig. 4,10 Body pitch Angular rate (q) vs Time (t) 
CBTT of Elliptical airframe; 2 Gees 
(0°,180°) 
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Fig. 4.11 Roll angular rate (p) vs Time (t) 

CBTT of Elliptical airframe; 2 Gees 
(0°, 180°) 
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Fig. 4.13 Roll angle { 4 ,) vs Time (t) 

CBTT of Elliptical airframe; 2 Gees (0®, 180°) 
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Fig. 4.14 



Pitch tail incidence (^p) vs Time (t) 

CBTT of Elliptical airframe; 2 Gees (0°/ 180°) 
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Fig. 4.15 Yaw tail incidence (6^) vs Time (t) 

CBTT of Elliptical airfframe; 2 Gees (0°/ 180 ) 
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Fig. 4.16 Roll tail incidence (5^^) vs Time (t) 

CBTT of Elliptical airframe; 2 Gees (0°/ 180 
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cross-coupling was removed. The major influence of the 
aerodynamic cross-coupling in the elliptical airframe 
response has been to decrease the overshoot in the achieved 
maneuver plane acceleration resulting from the second 
guidance command. 

In addition to aerodynamic cross-coupling the 
kinematic cross-coupling of -3p into a and inertial 
cross-coupling of pr into q were removed. Therefore, the 
only cross-coupling which exist in the airframe dynamic 
model are the kinematic coupling of ap into 6 and the 
inertial coupling of -qp into r. There is also the 
autopilot cross-coupling of the coordinating command from 
the Roll to Yaw channel. 

The equations, which represent the modified 
aerodynamic model by removing the referred cross-coupling 
are same with those as in paragraphs IV.B.l, IV. B. 2, IV.C.l 
and IV, C. 2, with the only difference in the following 
equations; 

1. Pitch Aerodynamic Model 



^ d Cm flyy C^€.g/s2.C^) 

z i Cm 



(IV. E. 1-1) 



and 



a z c)4- 

Z 9 - ( K ^ ^ Cfj /w') 

Z C| - O.C59^8 Cm (d^^/sec) 



4 



(IV. E. 1-2) 
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2 . Lateral Aerodynamic Model 

(Cn^.iS+Cnjy.^y) 

r-a04W52 + (dc^jsec^) 

(IV. E. 1-3) 



and 



P - 



5^-3 q,^cl 






Txx ' 

z 5^00^B5J^ (de^lsec^) 



(IV. E. 2-2) 



For purposes of analysis, a CSMP program (Appendix 
H) was written, using these equations. Figures 4.17 
through 4.24 show the time responses of the n„, n„ , a, e, 
q, r, S , 6„, fir,. The above figures and table II show the 

-3»/p»Y'R 

time constants and the overshoots of the variables. 

Figure 4.17 and Table II show that the achieved 
inertial maneuver acceleration satisfy the requirements for 
the time constants and overshoot, when the cross-couplings 
which were discussed before have been removed. 

Figure 4.25 (taken from [Ref. 2] shows the 
critical feedback paths which couple the pitch and yaw 
channels via missile roll rate. 

The aforementioned figures and plots verify the 
results of the [Ref. 2 ], 



F. CONTROL LAWS FOR CIRCULAR AIRFRAME 

The control laws used for the circular airframe are 

the same as those used in the linear studies, except for 
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Fiq. 4.17 Achieved Inertial Acceleration (n^) vs Time (t) 

CBTT of Elliptical airframe; Ideal 
Airframe dynamics; 2 Gees (0°,180°) 
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Fig. 4.18 Angle of Attack (a) vs Time (t) 

CBTT of Elliptical airframe; Ideal 
airframe dynamics; 2 Gees (0°, 180°) 
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Fig, 4.19 Sideslip angle ( 3 ) vs Time (t) 

CBTT of Elliptical airframe; Ideal 
airframe dynamics; 2 Gees (0°, 180°) 



145 





1 




010E6/SEC) 



O 




Fig. 4.20 Body pitch angular rate (q) vs Time (t) 
CBTT of Elliptical airframe;^Ideal^ 
airframe dynamics; 2 Gees (0 , 180 ) 
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Fig. 4.21 Body Yaw angular rate (r) vs Time (t) 
CBTT Elliptical airframe; Ideal 
airframe dynamics; 2 Gees (0°, 180®) 
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Fig. 4.22 Pitch tail incidence (fi ) vs Time (t) 
CBTT of Elliptical airframe Ideal 
airframe dynamics; 2 Gees (0°f 180°) 
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Pig. 4.23 Yaw tail incidence (5p) vs Time (t) 
CBTT of Elliptical airframe; Ideal 
airframe dynamics; 2 Gees (0°, 180°) 
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CBTT of Elliptical airframe;^Ideal^ 
airframe dynamics; 2 Gees (0 r 180 ) 
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Fig. 4.25 Critical Feedback loop for CBTT 
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TABLE II. CBTT PERFORMANCE OF ELLIPTICAL AIRFRAME 



VARIABLE VALUE TIME (sec) 



^Zrnin 


-2. 934 (gees) 


5.061 


^2max 


1 . 351 (gees) 


5.661 


"zmln 


-2. 934 (gees) 


5.062 


^Zmax 


0 . 937 (gees) 


5.053 


^m in 


1. 35 (deg) 


6.10 


a 

max 


6.75 (deg) 


3.10 


^m in 


2.22(deg) 


6.295 


^max 


6.75 (deg) 


3.10 


®min 


-0.9(deg) 


5.49 


®max 


1.49 (deg) 


6.05 


®min 


1. 53(deg) 


5.45 


®max 


1.07 (deg) 


5.91 


^min 


-10. 6 (deg/sec) 


5.639 


^max 


9. 5 (deg/sec) 


2.337 


% in 


-7 . 8 (deg/sec) 


5.328 


q 

^max 


9. 5 (deg/sec) 


2.337 


^m in 


-14. 115(deg/sec) 


6.600 


^max 


282.094 (deg/sec) 


5.387 


^min 


-14.29 (deg/sec) 


6.600 


^max 


201 . 928 (deg/sec) 


5.387 


^^min 


-2.32(deg/sec) 


5.90 


^max 


33 . 1 (deg/sec) 


5. 32 


^min 


-0 . 170 ( deg/sec ) 


7.036 


'^max 


37. 5 (deg/sec) 


5.29 



Note: 



CROSS- 

COUPLING 


T. <sec) 
Overshoot% 


T^Csec) 
Overshoo t% 


All 


0.48 


and 


0.41 and 


All 


0% 




11.25% 


ideal 


0.48 


and 


0.37 and 


ideal 


0% 




0% 



All 

All 

Ideal 

ideal 

All 

All 

ideal 

ideal 

All 

All 

ideal 

ideal 

All 

All 

ideal 

idea 1 

All 

All 

ideal 

ideal 



All the above referred values are 



T2 = 



62 percent time constant of achieved maneuver plane 
acceleration due to the first guidance command. 

63 percent time constant of achieved maneuver plane 
acceleration due to the second guidance command. 
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the gain (a) of the coordination branch of lateral control 
laws, as mentioned in the study for the elliptical 
airframe . 

The same 2 gees (0°, 180°) guidance commands and 
flight condition were used for direct comparison with the 
performance results of the elliptical airframe. 

However, to decrease simulation run time, the first 
guidance command and anti-gravity bias were applied at zero 
time with no missile or autopilot initial conditions. 

1 . Pitch Control Law 

The pitch control law for the circular airframe is 
the same as that use in the uncoupled autopilot shown in 
Figure 4.26. The command used in the autopilot is 
different than the one used in uncoupled channel and it is 
as follows: 



y\c- cose 



(IV. B. 1-1) 



as referred to in section IV.B.l 

The equations which represent the nonlinear pitch 

control law are the same as those in section III.C. of 

the linear pitch control law, using = -0.0387 

2 

(IV. F. 1-1) 
(IV. F. 1-2) 
(IV. F. 1-3) 



The equations are: 

i =-i50X+- 130 Hje 
%- Vh - 13531 (vic-x) 
Y, = -5V4 +6.573-(Hc-X) 



Sf^~ -O.m.Spt-O.3‘)058(V3+q/5T.3)-3.23O? (V3-qtaJ) 

(IV. F. 1-4) 
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Jp = -£8g. cfp + IOT95.3 *fpc (IV. F. 1-5) 



The state variables of the above equation are 
shown in Figure 4.26. 

2. Lateral Control Laws 

The lateral control laws for the circular airframe 
are the same as those used for the linear studies, except 
the modification in a^ (or a) which was referred to before. 



airframe. The equations which represent the lateral 
control laws are the same those of section III.C.l and 
III.C.2. and they are; 



The same command is used, as in the elliptical 



Yj : -5Vj + Li Ityg 




(IV. F. 2-1) 




(IV. F. 2-2) 



where : 




6 



and 




(IV. F. 2-3) 



r -8X1 -iTC 4> + iT.^ (iv.F.2-4 

Z -5 4 0.O5O33CX,-P/sT.3)-h0.?55 (Xf-P/sT-j) 



(IV. F. 2-4) 



(IV. F. 2-5) 
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X2 = -6X2 +0.0^2 ( P/573') 

-15^i2c +0.25 (yrX2)+<5(V^-X2') 
4 ^ 0^95.3 



(IV. F. 2-7) 



(IV. F. 2-6) 



(IV. F. 2-8) 



The above state variables are shown in Figure 



4.27. 

G. AERODYNAMIC MODELS FOR CIRCULAR AIRFRAME 

The aerodynamic models of Pitch and Roll-Yaw channels 
for the circular airframe are the same as those of the 
elliptical airframe and are shown in Figures 4.1, 4.2, the 
only differences which exist between them, are the 
aerodynamic mass parameters and nonlinear derivatives 
changes, as referred to in the (Appendices A and B) . 

Hence, the equations which represent the dynamic 
models are: 

1 . Pi tch Channe 1 



parameter for circular airframe, the above equations 
become ; 




(IV. C. 1-1) 



q ~-hL- + -53 : 1, 5 (Jc^/sec^) 



(IV. C. 1-2) 



a r C) + V<. (deg !sec) 



(IV. C. 1-3) . 



Substituting the values for aerodynamic and mass 



y\l^ - -2.05292 • Cf/ 
q r-.0.0<?^52 p.r+736.S60i|Cii( 



(IV. G. 1-1 ) 



(IV.G.1-2) 
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Actuator 




r- 

CN 

CX4 
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CBTT nonlinear lateral control law 



r 



if. 




a ~ Cj - 0.995H 0.0i^M52 p-^ (lv.G.l-3) 

Lateral Channels 

CCy^.’^ + -W / W (9€es) (iv.c.2-1) 

4 - K.viyg + a p/57.3 -r (devisee) (iv.c.2-2) 

(IV. C. 2-3) 

p r . € 4 ,S^') 

Ixx 

(IV. C. 2-4) 

<^>4 r p Cde0(5ec') (iv.c. 2 - 5 ) 



Substituting the values in the above equations, 
they yield: 

r 2.05292 (Cy^*^+ <^Y<Jy *^y) (IV. G. 2-1) 

4 r o.gs5^(Cy^.^4Cy^y-<^y) +0.0^5 p-a-r (iv.G. 2 - 2 ) 
Y r -6.0^7452 <^-P 4753.3837 (Cy4-^4 

(IV. G. 2-3) 



H. ANALYSIS OF NONLINEAR CBTT AUTOPILOTS FOR CIRCULAR 

AIRFRAME 

For purposes of analysis, a CSMP program was written 
(Appendix H) with the equations representing the 
aerodynamic models and control laws of the nonlinear CBTT 
autopilots for referred airframe. 

At first, in order to understand the effect of faster 
response, in the reduction of the effect of kinematic and 
inertial cross-coupling, the pitch acceleration error gain 
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0.0274 for slower 



Figure 4.26 was reduced to = 

response . 

The equation, which represents the acceleration 
compensator after the above reduction, is; 

s -(-|-+0 

Utilizing method of state-space representation with 
forcing term derivatives, the equation IV.H-1 yields; 

Y^-0.'35S0M•{->Ic-X) (iv.H- 2 ) 

Yi) = "5-Y + 4.t53 ■('Vlt'X) (lv.H-3) 

Using the above equations, instead of the equation 
IV. F. 1-2 and IV. F. 1-3, a new CSMP program (Appendix I) was 
written . 

Figure 4.28 shows that the achieved maneuver plane 
acceleration has a good response to the first command which 
is applied at zero time. The missile with the circular 
airframe moves upward like a skid-to-turn missile, as did 
the elliptical airframe missile, because the motion is in 
the desired maneuver direction and therefore the roll 
channel is not commanded. 

The time constant of the achieved maneuver response 
is equal to 0,49 seconds. The achieved maneuver plane 
response due to the second guidance command applied at 3 
seconds, shown in Figure 4.28 is reacting differently to 
the kinematic and inertial cross-coupling than the 
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Fig. 4,28 Achieved Inertial Acceleration (n^) vs Time (t) 

CBTT of Circular airframe; 2 Gees (0°, 180°) , 

K, = -0.0274 
^1 
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(iRf ■ ; 



elliptical airframe in Figure 4.5. Rather than overshoots 

and undershoots, a slowing transient starts at 3.5 seconds. 

Figure 4.29 shows that the overshoot in the body-fixed 

pitch acceleration, due to the kinematic and inertial 

cross-coupling during the second command, is 60.8 percent 

which is substantially more than it was for the elliptical 

airframe and occurs much sooner at 3.66 second. 

Figure 4.29 through 4.39 show the time responses of 

the ny , a, 5, q, p, r, 4> , «p, « and Above figures 

^c 

and table III show the characteristic of the responses of 
the variables for the circular airframe, when the = 

-0.0274. 

The referred results of this study verify the results 
of [Ref. 2] . 

I. EFFECT OF INCREASING PITCH CHANNEL SPEED OF RESPONSE 
In order to reduce the effect of kinematic and 

inertial cross-coupling during the second guidance command, 
the response of the Pitch channel of the CBTT autopilot for 
the circular airframe was made faster as shown in section 

II. This was accomplished by increasing 

the acceleration error gain K. = -0.0387. 

2 

Using the above gain K. the equations represent the 

2 

acceleration compensator are the IV. F. 1-2 and IV. F. 1-3. 

The effect of the change in on achieved body fixed 
acceleration (Figure 4.40) results in the achieved maneuver 
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0.0 0.5 ].0 1.5 2.0 2.5 5.0 3.5 ^.0 4.5 5.0 5.5 5.0 

TlflElSECJ 

Fig. 4.29 Achieved Body-fixed Acceleration (n„ ) vs 

Time (t) CBTT of Circular airframe; 2 Gees (0°, 
180°) , = -0.0274 
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Fig. 4.30 Achieved Body-fixed Acceleration {ny ) vs 

B 

Time (t); CBTT of Circular airframe; 2 Gees (0°, 
180°) , K, = -0.0274 
^1 
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Fig. 4.31 Angle of attack (a) vs Time (t) 
CBTT of Circular airframe; 

2 Gees (0°, 180°), K, = -0.0274 

^1 
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Fig. 4.32 Sideslip angle ( 0 ) vs Time (t) 
CBTT of Circular airframe; 

2 Gees (0°, 180°), K, = -0.0274 

^1 
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Fig. 4.33 Pitch Body angular rate (q) vs Time ( t) 
CBTT of Circular airframe; 2 Gees 
(0^^, 180°) , K, = -0.0274 
^1 
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Fig. 4.34 Yaw Body angular rate (r) vs Time (t) 
CBTT of Circular airframe; 

2 Gees (0°, 180°), = -0.0274 
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Fig. 4.35 Roll Body angular rate (P) vs Time (t) 
CBTT of Circular airframes; 

2 Gees (0°, 180°), K, = -0.0274 

'^l 
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Fig. 4,36 Roll angle (<j») vs Time (t) 

CBTT of Circular airframe; 

2 Gees (0°, 180°), K. = -0.0274 

^1 
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Fig. 4.37 Pitch tail incidence (6p) vs Time (t) 

CBTT of Circular airframes; 

2 Gees { 0 °, 180°), K, = -0.0274 

^1 
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Fig. 4.38 Yaw tail incidence (6^) vs Time ( t) 
CBTT Circular airframe; 

2 Gees (0°, 180°), =-0.0274 
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Fig. 4.39 Roll tail incidence ( 6j^) vs Time (t) 
CBTT Circular airframe; 

2 Gees (0°, 180°), K, =-0.0274 

^1 
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TABLE III 



CBTT PERFORMANCE OF CIRCULAR AIRFRAME 



VARIABLE 


; VALUE 


TIME(sec) 


CONDITION 


(sec) ; 
oOershoot% 


T2 ^ > 
overshoo t% 


“^Zmax 


1.092(gees) 


5.10 


^A 


0.46 


and 


0.52 and 


*^Zmin 


-3.044 (gees ) 


3.060 


% 




% 


^Zmax 


1. 109 (gees) 


4.540 


^A 

^2 


0. 34 


and 


0.62 and 


*^Zmin 


-3.125(gees) 


0.825 




5% 




0.05% 


^Zmax 


1 .002 (gees) 


5.040 


^A 

^3 


0.34 


and 


0.53 and 


*^Zmin 


-3. 125 (gees) 


0.825 




5% 




0% 


^max 


11.15(deg) 


1.432 


^A 








^min 


0 (deg) 


0 








^max 


11.489 (deg) 


0.862 


^A 

^2 








^ i 

min 


0.053(deg) 


3.656 










^max 


11.488 (deg) 


0.862 










^min 


0 ( deg) 


0 










®max 


4.25(deg) 


3.461 










®mln 


-0.399(deg) 


4.286 










®max 


3.633(deg) 


3.450 


•'a 

2 








®min 


-0.326(deg) 


4.290 










®max 


5.114 (deg) 


3.608 


^A 

3 








®min 


-0.414(deg) 


4.620 










^max 


21 . 39 (deg/sec) 


0.285 


^A 

^1 








^min 


-8 . 84 (deg/sec) 


3.220 








^max 


29.562 (deg/sec) 


0.261 










^mln 


-14 . 323 (deg/sec) 


3.225 










*^max 


29 . 562 (deg/sec) 


0.261 










^min 


-17. 575 (deg/sec) 


3.251 










^max 


278.5(deg/sec) 


3.394 










^raln 


-14.166 (deg/sec) 


4.6 










^max 


279.487 (deg/sec) 


3.4 










^m in 


-14.071 (deg/sec) 


4.6 










^max 


277 . 574 (deg/sec) 


3.4 


3 








^min 


-13.798 (deg/sec) 


4.597 










r 

max 


22 . 66 (deg/sec) 


3.28 


''a 

''l 








*^min 


-1.456 (deg/sec) 


4.56 








r 

max 


21.406 (deg/sec ) 


3.262 
















2 








^rain 


-1 . 322 (deg/sec) 


4.537 










^max 


23.159(deg/sec) 


3.287 


'•''3 








r . 
min 


-1.364 (deg/sec) 


4.818 
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Time (t) CBTT of Circular airframe with faster 
responding Pitch channel; K. = -0.0387, 2 Gees 
(0°, 180°) ^2 
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plane acceleration response (Figure 4.41). Figure 4.41 
shows that the achieved maneuver plane acceleration 
response during the first guidance command has improved. 

Table III shows that increasing the resulted in a 
decreasing in maximum sideslip angle 3.66 and a slight 
decrease in yaw angular rates (maximum 21.4 degrees). 

J. EFFECT OF INERTIAL AND KINEiMATIC CROSS-COUPLING IN 
PITCH CHANNEL 

The cross-couplings, (-0P) into a and (rP) into q , 
were removed in order to assess their effect on 
performance. Hence, the equations for a and q in the pitch 
aerodynamic model becomes: 

d - q -f K. Tlig (devisee) (iv.j-i) 

q liyy dv.j-2) 

Replacing the values to the parameters of the above 
equations, they yield: 

a = c|-o.q?SH Cfj (iv.j-3) 

q =m.wcn uv.j-4) 

where C^^, Cjj^ nonlinear functions which vary with a and 
5 referred to in (Appendix B) . The lateral aerodynamic 
cross-coupling was retained. For purpose of analysis, a 
eSMP program (Appendix I) was written. 

Comparing Figures 4.42 and 4.43 with figures 4.40 and 
4.41 it is concluded that the undersirable transients are 



gone 
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Fig. 4.41 Achieved body-fixed Acceleration (n„ ) vs Time 

(t) CBTT of Circular airframe; Faster responding 
Pitch channel; 2 Gees (0°, 180°) = -0.0387 

^^2 
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Fig. 4.42 Achieved Inertial Acceleration (n^) vs Time (t) 



CBTT of Circular airframe; Faster Pitch 
channel; inertial and kinematic cross-coupling 
into Pitch channel removed; 2 Gees (0°, 180°)/ 



K, = -0.0387 
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. 4.43 Achieved Body-fixed Acceleration (n„ ) vs Time 

(t); CBTT of Circular airframe; Faster Pitch 
channel; inertial and kinematic cross-coupling 
into Pitch channel removed; 2 Gees (0^, 180®) r 
K, = -0.0387 

^^3 
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Table III shows that the achieved maneuver plane 
acceleration time constant for the second guidance command 
has decreased. 

Figures 4.44 through 4.50 show the time responses for 
a, B, q ,r, 6p, <Sy, and Bp when the is equal -0.0387 

and the kinematic and inertial cross-coupling of the Pitch 
channel have removed. 

K. CONCLUSIONS 

1. The maneuver plane responses have transients, 
which may have to be reduced, are caused by inertial and 
kinematic coupling between pitch and yaw dynamics. 
Transients cause excessive overshoots and undershoots in 
achieved maneuver plane acceleration for elliptical 
airframe (Figure 4.6) and excessive slowdown in the speed 
of response for circular airframe (Figure 4.28, 4.40). 

2. The result of the nonlinear 3-D performance study 
verify the linear study of chapter III. 

3. The results of this study verify the results of 
[Ref. 2]. 

Further analysis and conclusion about a, q, r, p 
and 6 are discussed in [Ref. 2]. 



179 



mecccESi 




Fig. 4.44 Achieved Body-fixed Acceleration (ny ) vs Time 

(t) CBTT of Circular airframe; Faster pitch 
channel; inertial and kinematic cross-coupling 
into pitch channel removed; 2 Gees (0^, 180*^), 
K- = -0.0387 
^3 
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Fig , 




4.45 Angle of Attack (a) vs Time (t) 

CBTT of Circular airframe; Faster Pitch 
channel; inertial and kinematic cross-coupl* 
into Pitch channel removed; 2 Gees (0°, 180 



K- = -0.0387 
^3 
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Fig. 4,46 Sideslip angle ( g) vs Time (t) 

CBTT of Circular airframe; Faster Pitch 
channel; inertial and kinematic cross- 
coupling into Pitch channel removed; 

2 Gees (0°,180°), K, = -0.0387 

3 
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Fig . 




4.47 Pitch Body angular rate (q) vs Time (t) 

CBTT of Circular airframe; Faster Pitch 
channel; inertial and kinematic cross-coupl’ 
into Pitch channel removed; 2 Gees {0°, 180 
K, = -0.0387 
3 
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Fig . 




4.48 Yaw Body angular rate (r) vs Time (t) 

CBTT of Circular airframe; Faster Pitch 
channel; inertial and kinematic cross-coupl" 
into Pitch channel removed; 2 Gees (0°, 180 ; , 



K, = -0.0387 
^3 
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CBTT of Circular airframe; Faster Pitch 
channel) inertial and kinematic cross-coupling 
into Pitch channel removed; 2 Gees (0 / 180 ) , 



K, = -0.0387 
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CBTT of Circular airframe; Faster Pitch 
channel; inertial and kinematic cross-coupling 
into Pitch channel removed; 2 Gees (0°/ 180 °) , 
= -0.0387 
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V. 



MINIMIZATION OF THE KINEMATIC AND INERTIAL 



CROSS-COUPLING EFFECTS 



A. GENERAL 

The nonlinear CBTT autopilot for both airframes, 
consisted of the uncoupled autopilots discussed in section 
II, which coupled via inertial, kinematic and aerodynamic 
couplings . 

The effects of the cross-coupling were discussed for 
both airframes in section IV. 

In this section the kinematic, inertial and 
aerodynamic coupling effects will be sumarized. The 
kinematic and inertial effects in the pitch channel will be 
analyzed . 

The minimization of the above coupling and their 
effects will be studied, using feedbacks of angle of attack 
and rate of angle of attacks. 

B. CROSS-COUPLING OF CBTT AUTOPILOTS 

As discussed in previous sections, the four 
cross-couplings, between the autopilot channels for both 
airframes are: 

1 . Aerodynamic C ro ss-Coupling 

The aerodynamic cross-coupling is the coupling 
between the Roll and Yaw nonlinear channels, via the 
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C„ . Co and C (Figure 4.2). The major effect of the 
6 

Y R 

above coupling in the elliptical airframe response is to 
decrease the overshoot in the achieved maneuver plane 
acceleration resulting from the second guidance command, 
(i.e 2 gee (180°), the other effects on the missile 
variables were small but in a direction which improves 
performance (i.e decreased missile body angular rates, less 
control surface motion, less sideslip variations). In the 
circular airframe, the effects of the Cj^ , C are 

negligible, the critical coupling coefficient is the C, 

which effects the stability in the coordination branch, 

which can be improved by reducing the effect of C, on 

5 

Y 

missile roll angular acceleration. (equation IV. C. 2-4) 

2. Cross-Coupling Between Roll-Yaw Autopilots for 
Coordination Motion 

This cross-coupling branch, shown in Figure 4.4, 
has been added between the yaw and roll control laws, in 
order to provide coordinated motion between these two 
channels. For purposes of stability the coordinating 
branch gain K^p was selected for the elliptical airframe 
Kyp = 1 and for the circular airframe Kyp = 0.458. The 
elliptical airframe is more stable than the circular with 
the above selected Kyp. 
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3. Kinematic and Inertial Cross-Coupling 



The kinematic and inertial cross-coupling exist 
between the pitch and lateral nonlinear aerodynamic models 
of the CBTT autopilots. There exist the following 
cross-coupling paths between pitch and lateral channel. 

a. Kinematic cross-coupling path of (-^p) into a 
in the pitch aerodynamic model (Figure 4.3). 

b. Kinematic cross-coupling path of (ap) into 6 
in the lateral aerodynamic model (Figure 4.4). 

c. Inertial cross-coupling path of (pr) into q in 
the lateral dynamic model (Figure 4.4). 

d. Inertial cross-coupling path of(-qp) into r in 
the lateral dynamic model (Figure 4.4). 

The efffects of the inertial and kinematic 
coupling are the transients in the achieved maneuver plane 
acceleration response, and have the form of overshoots and 
undershoots discussed in section III, IV. Transients may 
also result in a slower maneuver plane acceleration 
response and are less pronounced at higher acceleration 
levels . 

Coupling transients may be reduced by increasing 
pitch stability, by decreasing maximum roll rate (section 
IV) , by improving the autopilot coordination technique to 
minimize sideslip rate S (section IV) , by increasing the 
effects of stabilizing lateral aerodynamic coupling 
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{negative C ) r by increasing pitch channel speed of 



response (section IV) and by increasing the pitch stability 

C„ with feedbacks of both angle-of-attack and rate 
m ^ 

a 

of angle-of-attack for the pitch autopilot. 

C. KINEMATIC AND INERTIAL CROSS-COUPLING IN PITCH 
AUTOPILOT 

The kinematic and inertial cross-couplings in the 
pitch autopilot, are via the pitch angular rate q and the 
rate of angle-of-attack (a) (Figure 4.1)). 

The equations for q and a (figure 27) are; 



Removing Cj^^ from the parenthesis in (V.C.l) and 
"*a 

combining with (V.C.2) to eliminate q, the (V.C.l) results 




(V.C.l) 



(V.C.2) 



where : 



rP/57.3 = Inertia Cross-Coupling 

-3P/57.3 = Kinematic Cross-Coupling 



in; 




(V.C.3) 
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Let ; 



and 









Ivy 



>na 



-i^cc. 



Taking the Laplace transform of (V.C.3) and solving 
for a, the equation becomes: 

-A f (rP/sT.S)- S-(^) r\ 

Ct= r • f»f2*<3P I (V.C.4) 

H\l I Ki j 

Kl 

Increasing the magnitude of will decrease both K 2 

‘a 

and 1/K^ resulting in the reduction of cross-coupling 
effects on the angle-of-attack ( o) and the influence of 
in the a . 

In the section III the equation of n was derived: 

r (-q,^/w) ‘ ) (V.C.5) 

which involve (a). Therefore^ the reduction of inertial 

and kinematic effects on (a) , results in the reduction of 

the effects in the achieved body-fixed acceleration (n_ ) . 

B 

In the nonlinear studies the equation for the achieved 
inertial acceleration was derived as: 

-COi4>+'YJY^*Sl'n4) (V.C.6) 

It is apparent from the equations (V.C.4), (V.C.5) and 

(V.C.6) , that an increasing of C^^^ results in the reduction 

"“a 

of kinematic and inertial effects (i.e transients in the 
response) in the achieved inertial acceleration (n 2 ) » which 
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it is desired in order to meet the requirement of time 
constants (Tj^, ^ 2 ) and percents of overshoots and 
undershoot, as referred in Appendix C. 

D. MINIMIZATION OF THE KINEMATIC AND INERTIAL 

CROSS-COUPLING EFFECTS IN ELLIPTICAL AIRFRAME 

As discussed in the nonlinear studies, using the pitch 
angular rate q as feedback and under the influence of the 
inertial and kinematic effects, the achieved inertial 
acceleration (n 2 ) for both airframes doesn't meet the 
requirements of the time constants, overshoot and 
undershoot, denoted for the linear autopilots, referred to 
in Appendix C. 

Figure 4.5 and Table II show that the achieved 
inertial acceleration (n^) for elliptical airframe with 
feedback the pitch angular rate (q) has overshoot 11.25% 
and = 0.48 sec, ^2 = 0.41 sec. 

Figures 4.40 and Table III illustrate that the 
achieved inertial acceleration (n 2 ) for circular airframe, 
with feedback the (q) has overshoot 5% and r=-.46 sec, 

■^2 = -.62 sec. Hence the (n 2 ) for both airframes is not 
acceptable. 

In order to minimize the inertial and kinematic 
cross-coupling effect, namely to have smoother transient 
response of n 2 which will meet the requirements of time 
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constants and overshoot, the magnitude of was increased, 

a 

utilizing as feedback the angle-of-attack and its 
derivative . 

In order to increase the magnitude of C„ the 

^ in 

a 

stability augmentation theory was used. 

Considered as stability derivative augmentation ®— ►'Sp 
and increase the magnitude of and directly, as 

referred in [Ref. 6], the increases in the Ma is 

9 

essentially an increase in and hence the inertial and 

a 

kinematic effects on the angle-of-attack are reduced. 

In the pitch control law (Figure 3.1) a SAS exists 
using as feedback the pitch angular rate (q) . In order to 
design the SAS with the angle of attack (a) as feedback, 
the root locus and lead-lag compensation techniques were 
used . 

1 . Stability Augmentation System (SAS) With Feedback 
of the Pitch Angular Rate (q) . 

In Figure 5.1, the stability augmentation system 
is shown, with feedback of pitch angular rate q where the 
transfer function of the q/<S_ for the uncoupled pitch 
dynamic model, given in section II. C. The closed-loop 
transfer function of the above SAS is; 



k(T+0^ST3) (-MJ.35 S- I.TflS) 

; ; (v.D.i-i) 

S (5?3) (:^ + 1) ( S^+0.-(8JS-ii.2T9) 
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(V.D.1-2) 



Manipulating equation V.D.1-1, it results; 



6 ( 5 )= 



0.0053* 1.0009 -sV 0.(ll*i^- 1J18 • S + K (-5.f69 S- tfj. 321S-1>^11S) 



The poles and zeros of the SAS are as follows; 

a. = -188.43 

b. ?2 = -3.45 

c. P3 = 3.269 

d. = 0.0 
and 

6 . Z3 = —8 . 0 

f. Z2 = -0.186 

The other two zeros are at infinity. 

The root locus of the system is shown in Figure 

(5.2). 

The K was selecting at the value of -3.07, because 
at this value, the system has two dominant roots at P^'2 ~ 
-9 . 212 ± j7. 606 which dominate the response of the system, 
having the following characteristics; 

a. 0 = 

b. J z COS 9 = 0.733 

c. W- Wrif =1615 (roidlscc) 

d. UUt) = 

where ; 

; : damping constant 

00 ; frequency of system 

u)^; natural undamped frequency 
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2. SAS with Feedback of the Angle of Attack (a) and 
cate of angle of attack (a) 

The SAS with feedback the pitch angular cate q 
(Figure 5.1) and K = -3.07, is stable and the dominant 
roots gives damping ratio 5 = 0.73 in the response of the 
system . 



To minimize the coupling effects the angle of 
attack (a) was selected as feedback in the previous SAS. 

It is desired the SAS to have the same behavior as before, 
for purposes of stability and response of whole the pitch 
channel . 



Pole zero cancellation technique is used, to 
obtain the same poles and zeros, as in the case with 
feedback the pitch angular rate q. 

For above purposes, a compensator was used in the 
SAS, as shown in Figure 5.3. 

The transfer function of the compensator is: 



~ -0.02 S - q(.35^ (V.D.2-1) 

The closed loop transfer function of the SAS is as 



follows : 








Yx(4-*hY(' «.3 y 


/-0.02S -i^f.35T \ 


\-ao2-s-m.35 ? ) 




\5^+OJ8fS*iJ.2?9 j 



H 




\( 5T.5 


U-O02.S-i(i.35T \ 






jVS*+0.<81-5- 11.279 ) 
(V.D.2-2) 


where the 


transfer function a _ 


-0 w s - 


is 


derived in 


6 

p 

the section II. C. 3. 
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Pitch channel; CBTT of Elliptical airframe 



Rearranging and manipulating the equation 
(V.D.2-2), it becomes: 

_ U (5.9131 sV i229.38l -5 4- i82^S.O^) 

' -O.OOOl-S®+ a22l5 siiH.399 5^-</.7?9 5^- 466.466 5 + K (al034 S^.2I<J.6 sV<miS*3fJ.9) 

(V.D.2-3) 

The poles and zeros of the system are: 

= -1877.931 
p 2 = -207.391 
P 3 = -3.447 
P 4 = -3.271 

P3 = 0.0 

and 

Z 3 = -2067.838 

Z 2 = -8.0 

Z 3 = -0.186 

The other two zeros are at infinity. 

The root locus of the system is shown in figure 

5.4. 

Selecting the K = -3.07, the system has the same 
dominant roots P ^^^2 = -8.249 ±j7. 577, as in section V.D.l. 

3. Acceleration Loop 

Considering the achieved acceleration loop in the 
Pitch channel, the modification of the SAS, effects the 
transmission gain of the above loop. 
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In order to compensate the effect in the 



transmission path, a compensator was added as shown in 
Figure 5.5. 



inverse 


The transfer function of the compensator is the 

of the V.D.2-1: 

^ . 1 -0.02•S-l^i35? 

^ Gcs)~ (V.D.3-1) 

Modifying in this way, the transmission path of 



the acceleration loop, the gain of path, is the same as in 
the case of par. V.D.l. 



4. 


Measuring of the Minimization in the Kinematic and 

Inertial Cross-Coupling Effect 

In order to measure the minimization a CSMP 


program 


(Appendix J) was written using the equations of the 



nonlinear aerodynamic models, the equations of nonlinear 



lateral 


control law and the equations of the modified pitch 


control 


law. 

The equations of the modified pitch control law 


are the 


following : 

a. Acceleration lag network 

X= -150X4150 (II. C. 4-2) 

b. Acceleration second order compensator (Figure 


5.5) 


V _ (-0.00s).(-0.03*s-m.35?') 

(ur»c-x') (4.40(-‘)i35f5-r.wO 
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utilizing inverse Laplace transformation and 



rearranging the above equation, one has: 

Y+4H8S Y = -0.000232 (lJfi1c-x)-0.W6*(i iJVX) (V.D.4-2) 
Using space representation of a system, in which 
the forcing function involves derivative 
[Ref. 5, pp 675-678], one obtains: 

- V'(2-0.00025?5-\ +0.000232 -X (V.D.4-3) 

W 2 Z-J.118??W|-6.i86+V«^-0.53J2 \+0.^im X (V.D.4-4) 

V Z (V.D.4-5) 



c. Compensator of Actuator 
<^Pc . 

(y-a|si3)" s (-0.02-5-m.357) 

Using the method for state representation and 



(V.D.4-5a) 



rearranging the equation V.D.4-5 yields: 

= ip2 +K- 532.332X<0^-(Y-al5T.5) 



(V.D.4-6) 



-206+35 <JP2- X. -H00.926Xi0^(N-a|5+3) (V.D.4-7) 

5[>^z (Jp, - K. 2S8.i^8(Y-a|5+5) (v.D.4-8) 

d. Equation of Actuator 



-I8g.M-^p+ H0795.3-ipc (iv.B.l-4) 

Selecting K=3.07, for the compensator of SAS, with 
feedback (a) the time constant (i.e time constant of 
achieved maneuver acceleration due to the first command) is 
0.45 sec, the time constant (for second command) is 0.41 
sec and the undershoot is 20.73 percent as shown in figure 



5.6. 
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Modified Pitch channel; CBTT of Elliptical 
airframe; 2° (0°, 180°), K = -3.07 



The above response of the inertial achieved 
maneuver acceleration does not have, smooth transient, 
because of the large undershoot. 

Due to the above reason, another K was selected 
equal to -19.37, for minimization of the coupling effects. 

With this selected K the system has the following 
poles, shown in figure 5.4 



P^,2 = -98±j96.239 
= -1879.06 
= -8.596 
= -0.189 

With the same zeros as section V.D.2, the 



dominant poles of the system are 
characteristics of the response are; 



a. 9 = 44.24 

b. C = 0.716 

c. (0 = 93.775 rad/sec 

d. oj^= 134.412 rad/ sec 

It is shown in known (Ref. 6) that the can be 



expressed as: 




(V.D.4-9) 
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and 






(V.D.4-10) 



where : 

y: specific weight of air 1.4 

I: moment of inertia 

P: Local pressure 

M: Mach number 

A; Area 

Comparing the co^ = 134.412 for the K = -19.37, it 

is apparent that it is much larger than the = 11.203 for 

K = -3.07 and hence according to the equations V.D.4-9 and 

V.D.4-10, the increases essentially, resulting in 

™a 

minimization of kinematic and inertial cross coupling 
effects. 

Figure 5.7 shows, that with K = -19.37 and 
feedback in SAS the (a) , the achieved inertial maneuver 
acceleration has = -.49 sec, ^2 ~ “‘38 sec and overshoot 
about 0%. 

The above analysis illustrates that the selection 

of K = -19.37, instead of K = -3.07, keeps the damping 

ratio about constant c = 0.71 (C = 0.73 for K =-3.07) and 

increases the w , which essentially increases the , 

n' ■* m 

a 

resulting in very smooth response of the achieved maneuver 
acceleration . 



206 




\ 






LEGEND 



□ 




207 



airframe; 2^ (0^,180^), K = -19.37 



Using the same particular compensators in the 
Pitch channel shown in Figure 5.5 and feedback the rate of 
angle of attack (a ) , the Pitch channel becomes 
uncontrolable for both selection of K = -3.07 and K = 
-19.37, because the pitch tail incidence exceeds the limits 
of ±10°. 

Also, selecting the K = -19.37 in SAS of Figure 
5.1 with feedback the pitch angular rate q the pitch 
channel becomes uncontrolable again, due to the exceeding 
of the limits in 

Therefore, for minimization of the kinematic and 

inertial effects (i.e improving of the transient response 

of the achieved maneuver acceleration) for this particular 

airframe and flight conditions, the Cjjj was increased, 

a 

using the aforementioned compensators in the pitch control 
law (Figure 5.5) and feedback the angle-of-attack (a). 

Comparing the achieved inertial acceleration (n 2 ) 
with feedback (q) in the Pitch channel (Figure 4.41) and 
the achieved inertial acceleration (n 2 ) with feedback (a) 
and K= -19.37 (Figure 5.7), it is concluded that for these 
particular airframes and flight conditions the feedback of 
angle-of-attack (a) is more desirable, than the feedback of 
(q) and (a) , because in the achieved inertial acceleration 
has been reduced the kinematic and inertial effects, 
meeting with this way the requirement of the time constants 

and overshoot, denoted in the linear studies. 
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E. MINIMIZATION OF KINEMATIC AND INERTIAL COUPLING IN 

ELLIPTICAL AIRFRAME 

Making a comparison in Figure 5.8, which shows the 
kinematic cross-coupling (-p5) into a, in the pitch channel 
with feedback pitch angular rate (q) and Figure 5.10 which 
shows the -pS in the modified pitch channel (Figure 5-5) 
with feedback the angle-of-attack , it is concluded (Table 
IV) that the kinematic coupling is increased from 
(-0°. 988-5°. 877) to (-2° . 124-7° . 383) , when the kinematic 
coupling effects are reduced by the modified pitch channel. 

Figure 5.9 shows the inertial coupling (pr) into q in 
the nonlinear Pitch channel, with feedback the pitch 
angular rate (q) . 

Figure 5.11 shows the inertial coupling for the 
modified Pitch autopilot with feedback the angle of attack 
(a) (Figure 5.5). 

Comparing the above figures, it is concluded (Table 
IV) that the inertial coupling is increased from 
(-1°. 692-168°. 317) to (-20°. 373 - 181°. 415), when the 
inertial coupling effects are reduced by the modified pitch 
channel discussed in previous section. 

Hence, the minimization of the kinematic and inertial 
effect does not mean necessarily minimization of the 
cross-couplings themselves. 
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Inertial Coupling (P.r) vs Time ^t) 
CBTT of Elliptical airframe; 2(0 r 18® 
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Kinematic Coupling (-P.3) vs Time (t) 

CBTT of Elliptical airframe; Modified 
Pitch channel for minimization of coupling 
effect 
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CBTT of Elliptical airframe; 

Modified Pitch channel for minimization 
of the coupling effect 



For minimization of the kinematic and inertial cross 



coupling, the same modified nonlinear pitch control law was 
used as that of section IV. D (Figure 5.5), except the 
acceleration compensator which was removed from the 
transmission path of the n^, as shown in Figure 5.12. 

For purposes of analysis a CSMP program (Appendix J) 
was written, using the following equation for the nonlinear 
pitch control law and the same equation for lateral control 
law and dynamic models, as in section IV. B and IV. C. 



where K = 19.37. 

Figure 5.13 shows the minimized kinematic coupling and 
Figure 5.14 shows the minimized inertial coupling. Figures 
5.13, 5.14 and table IV illustrate that the minimization of 
the kinematic and inertial coupling is 264% and 195.26% 
respectively from its values in the nonlinear CBTT 
autopilot for elliptical airframe. Figure 5.15 shows the 
achieved inertial acceleration (n 2 ) , with the Pitch channel 
modified for minimization of kinematic and inertial 
coupling (Figure 5.12). From this figure it is concluded 



X r - 150X + 

^ 0.5305 'Vic + 

+ v<. 532.382x^0^- (V-a|5?.3) 

r -2 o61[.9S.(5’p 2 - k- ^^oo.mx^o^(Y-al57•3) 

5?^ ~K. 25?.M3 (Y-al5?.3) 



(IV. B. 1-1) 



(IV.B. 1-2) 



(V.D.4-6) 



(V.D. 4-7) 



(V.D.4-8) 




(IV.B. 1-4) 
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Kinematic coupling (-P.0) vs Time (t) 
Modified Pitch channel for minimzation 
of kinematic coupling; CBTT of Elliptical 
airframe; 2 (0°, 180®) 
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Modified Pitch channel for minimization 
of inertial coupling; CBTT of Elliptical 
airframe 



TABLE IV 



KINEMATIC AND INERTIAL CROSS-COUPLING FOR 
ELLIPTICAL AIRFRAME 

CONDITION INERTIAL (deg/sec^) KINEMATIC (deg/sec^) 

MIN MAX MIN MAX 

Elliptical airframe 



with cross- 
coupling 
(section IV. C) 


-1.692 168.847 -0.9881 5.877 


Elliptical ideal 
airframe (section 
IV. E) 


-1.698 167.37 -2.124 7.383 


Elliptical airframe 
with modified Pitch 
autopilot 
K = -19.37 
(section V. 1) 


-20.373 181.415 -1.88 8.644 


Modified Pitch 
autopilot K = -30.4 
(section V.D. ) 


-7.263 169.371 -1.956 8.711 


Modified Pitch 
autopilot for 
minimization of 
cross couplings 
(section V.E) 

K = -19.37 


-0.0577 57.5322 -0.523 1.363 



Note: K is the gain of the added compensator in SAS of the 
Pitch autopilot. 
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Modified Pitch channel for minimizatin of 
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that the kinematic and inertial effects (i.e transient in 
achieved acceleration response) is increased. 

F. MINIMIZATION OF THE KINEMATIC AND INERTIAL 
CROSS-COUPLING EFFECTS IN CIRCULAR AIRFRAME 
As shown in Figures 4.28, 4.40 and 4.42 the achieved 
maneuver acceleration in the circular airframe has very 
smooth response, in regards to the overshoot or undershoot 
on it . 

The effect of the inertial and kinematic 
cross-coupling in the nonlinear aerodynamic model is that 
the response of the inertial achieved acceleration has a 
slowing transient starts at 3.5 seconds, as discussed in 
section IV. G. Hence, minimization of the effects means in 
the circular airframe improvement of time constant during 
the second guidance command. 

1 . Stability Augmentation System (SAS) with Feedback 
the q. 

Using the same method for minimization, as in 
section V.D, the SAS of the pitch channel for the circular 
airframe with feedback q is, as shown in Figure 5.16. The 
transfer function of the q/^p is given in section II. C. 

The closed-loop transfer function of the above SAS 
is ; 
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(V.F.1-1) 



(4s+ ( - 59.1088S - 6.<3«H5) 

2(S) = ^ 

K. (cfe+1^ (sT-3) (-5eaoS8 S 

(olj3 ■*'^) (s^+OJ^28 5+i^i^.33^6) 

Rearranging equation V.F.1-1, one obtains: 

K {20.i9-S^-i- 30 f 00-25 S^+Q038Qlf.6S-h^3aQ.?) 

' a03?-5'^+7.00‘|5^ + 3.GT9'$^ + 3<0.3<3-S +^^33 +K- (-<3.02<»-s^- 

(V.F.1-2) 

The poles and zeros of the SAS are: (1) 

P^= -188.438; (2) 92^3 = -0 . 07 39 ± j6 . 657 ; (3) = -0.143; 

(4) Z 3 = -6.55 and (5) Z 2 = -0.118. The other two zeros 
are at infinity. The root locus of the SAS is shown in 
Figure 5.17. The selected value of K is -15.6 and the 
dominant roots are P]^ 2 “ -10. 736±j7. 536 (Figure 5.17). 

With those dominant roots, the response of the system has © 
equals to 35.06 degrees; C equals to 0.818; (o equals to 
7.545 rad/sec and equals to 13.125 rad/sec. 

2. SAS with Feedback of the a 

For minimization of inertial and kinematic effect 
in the circular airframe, the same technique was used (i.e 
pole-zero cancellation technique) as in section V.D.2 and 
the resulted pitch control law with feedback of 
angle-of-attack (a), is shown in Figure 5.18. 

The transfer function of the added compensator in 
the SAS is given by equation V.D.2-2. 
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(V.D.2-2) 
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Fig. 5.17 Root locus of SAS with q/6 : CBTT of 
Circular airframe 
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The closed loop transfer function of the referred 



before system is given by equation V.D.2-3: 



^2= 





l-l 


(-59.1-5-6.‘39 \ 




/-0.039-S-59.1 \ 


J 


(-0.039 5- 59.1 J ' 


life*' ) ' 


S^+0.1*185*M9.33 / 



K- 5T.5- (-59,1 S- 6.99) (-0.039 5- S9.0 

(-^+l) (-0.039-5-59.i)(^-H) (s+0'<*<8 >S+«JI|.33) 



(V.D.2-3) 



where the transfer function a/ is given in the 
section II. C. 3, rearranging and minor manipulating the 
above equation, one obtains equation V.D.2-4: 



(20.38- 30?00.2S 5*+ Z0382 ¥. 65 4- 23642. 8 ) 

' -aoOf5-2.^f5-^f‘/.3^5-229.?5^f83</9-5-2630.^+tf(0.365'+535.«S4555»6+i(<3) 

(V.D.2-4) 

The poles and zeros of the modified SAS are, 

= -1499.623, ?2 = -188.53, P 3 , 4 = -0 . 07 39 ±j6 . 656 , P 5 = 
-0.142, = -1499.6, %2 = -6.55, Z 3 = -0.1182 and all 

other zeros at infinity. 

The root locus of the system is shown in figure 
5.19, it is apparent from above Figure, that for k=-15.6 
the dominant roots are the same as those in section V.F.l. 

Considering the transmission of acceleration n 

c 

path, as in the elliptical airframe, a compensator was 
added in the path, for compensation of the effect in it, 
which is created by the added compensator in the SAS, as 
shown in Figure 5.20. 
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Root locus of SAS with a/'5_ 
CBTT of Circular airframe " 
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3. Measuring of the Minimization in the Kinematic and 



Inertial Cross-Coupling Effect 

A CSMP program (Appendix K) was written using 
equations of nonlinear lateral control law (section 
IV. G. 2), of the aerodynamic models (section IV. 4) and the 
following equations of the modified nonlinear pitch control 
law; 

a. Acceleration lag Compensator 

X r -150X +^50 Ylig (II. C. 4) 

b. Acceleration Second-Order Compensator (Figure 
3 - 20 ) 

\ (y p ^ 3-1) 

5 (|-+^)• (-59.10S-S- 6.38^5) 

utilizing inverse Laplace transformation and 
rearranging the above equation, yields 

% + 5.1183 Y + 0. 591 y . 0.022 SZ ^ 5 

(V.F.3-2) 

Using state-space representation referred to in 
(Ref. 5; pp 675-678), one obtains: 

% - 0.022^2 (V.F.3-3) 

V\l 3 - Ka- 3M.?50 (v.F.3-4) 

%z -0.59U-W2-5.11g-\^3^.\(;(^.ia39(vx) (v.f. 3-5) 
^3“ NVi (V.F.3-6) 
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c. Compensator of Actuator in SAS 

S?c _ (V.F,3-7) 

(V ccl5T3) (•^+H)*(-0.0394-5-59.i(H') 

Utilizing, the aforementioned method for state 



representation and rearranging the equation V.F.3-7, it 
yields : 

- <^p2 4 K (^3896.0^3S)* (Y 3 - C(/5f3) ( V . F . 3-8 ) 

S?2 - - 1500 <5p2-y.3.3X1o\cg5t3)(V. F. 3-9 ) 
c5'p^-K.(52.7i|4) (Ya- a/ 573 ) (v.F.3-10) 

d. Equation of Autuator 



ijp = -)«84 -Sp+ 10T95.3-?(> 



{IV. G. 3-8) 



Selecting the K = -15.6 and Ka = -0.0387 

2 

(Figure 4.26) and feedback the pitch angular rate (q) , the 
achieved inertial acceleration n^ , has = -.34 sec, = 
0.62 sec and overshoot 5% (Figure 4.40, Table III). 

This response of the n^, is unacceptable 
because it doesn't meet the requirement on the time 
constant, referred to in Appendix C. 

Selecting K = -15.6 and = -0.0387 (figure 
5.20) and feedback the angle of attack (a) in the pitch 
channel the achieved inertial acceleration (Figure 5.21) 
has = -.35 sec, x = 0.6 sec and overshoot 7%, which is 
again unacceptable. 
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CBTT of Circular airframe; K = -15. 
= -0.0387 



(V.G.3-9) 



It is known (Ref. 6) that: 

— = V 



J_r \ 

' \| 



Hence, in order to reduce the it is needed 



to increase the C, 



m 



Selecting K = -70.88 and = -0.0387 



(Figure 5.20) with feedback the angle of attack (a) , the 
dominant roots of the SAS (Figure 5.13) are ~ 

-90 . 55 ±j87 . 95 and the characteristics of the system 
response are e = 44.16 degrees, c = 0.79 , u3 = 87.932 
rad/sec and = 126.212 rad/sec. With that selected k the 



z decreases from 0.818 (when K= -15.6) to 0.72 and the “ 



n 



increases from 13.125 rad/sec (when K = -15.6) to 126.212 
rad/ sec . 



Figure 5.22 shows, the achieved inertial 

acceleration n„ , when K = -70.88, K, = -0.0387 and feedback 
^ ^2 

the (a). The = 0.33 sec, = 0.48 sec and overshoot 
14.33 %. It is apparent that the ^2 decreases, which is 
desirable but the overshoot increases, because the damping 
ratio c decreases. 

In order to reduce the overshoot and to keep 
the time constants below the 0.5 sec, .the reduces to 
-0.0274 (Figure 5.20) and the gain K^p (Figure 4.27) of the 
coordination branch in the nonlinear lateral control law 
increases from 0.458 to 1.0. 

Using K^p = 1.0 (Figure 4.27) in the 
coordination branch in the lateral control law, K, 
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CBTT of Circular airframe; K = -70.88, 
K,, = -0.0387 



LEGEND 
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CBTT of Circular airframe; K = -15 
= -0.0274 



LEGEND 
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CBTT of Circular airframe; K = -70.88 
= -0.0274 



a 




LECENO 
NZ □ 



Time 



(t) 



CBTT of Circular airframe; K = -70,88, 
-0.0274, Kyp = 1.0 
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= -0.0274 in the acceleration compensator and K = -70.88 in 
the added compensator in the SAS with feedback the 
angle-of-attack (a) (Figure 5.20), the obtained achieved 
maneuver acceleration has = -.43 sec, '^2 ~ ”*37 sec 
overshoot 8.9%, as shown in figure 5.25 and table V. 

Using in the pitch channel with the 
aforementioned particular compensators, as feedbback the 
rate of angle of attack (a) and making all the above 
selection of K and K^p, the whole system becomes 

uncontrolable because the required tail incidence angle 
exceeds the limits of ±10°. 

Also using the selections of K, Kyp and 

i feedback pitch angular rate q (Figures 4.26, 4.27), again 
the whole system becomes uncontrolable, because the 
required tail incidence angle exceeds the limit of ±10°. 

The above analysis illustrates that the 
angle-of-attack is more desirable feedback than the (q) and 
(a) , in the nonlinear pitch control law for these 
particular circular airframe and flight conditions, because 
it reduces the kinematic and inertial coupling effect, such 
that the achieved inertial acceleration (n 2 ) meets the 
requirements of the time constants and overshoot, referred 
to in the linear studies. 
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TABLE V 



CHARACTERISTIC OF THE INERTIAL ACHIEVED ACCELERATION ( 02 ) 
FOR CIRCULAR AIRFRAME 



CONDITION 






(sec) 


T 2 (sec) 


overshoot 


^A = 


0.0387, 


K = 


-15.6 


0.35 


0. 6 


7 


> 

It 


0.0387, 


K = 


-70.88 


0.33 


0.48 


14.33 


II 

< 


0.0274, 


K = 


-15.6 


0.435 


0.51 


8.6 


> 

It 


0. 0274, 


K = 


-70.88 








II 


1 






0.43 


0.37 


8.9 


Note 


; K,; pi 


tch 


autopilot 


acceleration gain in 


acceleration 



compensator 

63 percent time constant of achieved maneuver 
plane acceleration due to first command 
63 percent time constant of achieved maneuver 
plane acceleration due to second command 
K„p: gain of the coordination branch 

K; gain of the added compensator is the SAS of pitch 
autopilot 



1 * 



2 * 
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G. MINIMIZATION OF KINEMATIC AND INERTIAL COUPLING IN 

CIRCULAR AIRFRAME 

As in the elliptical airframe comparing the values of 
the inertial and kinematic cross-coupling for the circular 
airframe, namely, the (P.r) and (-p®) , referred to in table 
VI, it is concluded for this particular analysis that as 
the coupling effects in the nonlinear Pitch autopilot are 
reduced, the cross-couplings themselves are increased. 

In order to minimize the cross-couplings, the pitch 
control law was modified, as shown in Figure 5.26. 

For purposes of analysis, a CSMP program (Appendix K) 
was written, using the following equations for the 
nonlinear modified pitch control law, lateral control law 
and dynamic models, as in sections IV. G and IV. H. 

X - -150X + 150 > (II. C. 4) 

2. r W 2 - ' 0.022S8(>|(.”X^ (V.F.3-3) 

3. W 2 r Wa- K a- 314.850 (tIc-x) (v.f.3-4) 

W3--0.591W2'5.H8TV\/3+Ka-1?3.385{vx) (v.f.4-5) 

5. Y/i (V.F.4-6) 

6. -O.Hq3-Sp^-0.3i|058(Vot/5T*3)-2230^-(y3-cx/5?.3) 

(V.F.4-7) 

which is the equation IV. G. 1-4, except of the 
feedback, which is the angle of attack (a) instead of the 
pitch angular rate (q) . 

7. r + ^0?<^5.3•6pc (IV. G. 1-5) 
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TABLE VI 



KINEMATIC AND 


INERTIAL COUPLING FOR CIRCULAR AIRFRAME 


CONDTION 


INERTIAL (deg/sec^) KINEMATIC (deg/sec^) 



MIN MAX MIN MAX 



Circular airframe 
with cross- 
coupl ings 
(section IV. G) 


-a. 008 92.794 -1.132 0.012 


with cross- 
coupling and 
faster response 
K = -0.0387, 

K = -15.6 
(section V.F) 


-0.0768 100.573 -2.04 0.0353 


with cross- 
coupling and 
faster response 
K. = -.0387, 

K = -70.88 
(section V.F) 


-0.0212 103.553 -1.39 0.0436 


with cross- 
coupling and 
slower response 
K = -0.0274, 

= -15.6 
(section V.F) 


-0.107 118.294 -2.97 0.0387 


with cross- 
coupling and 
slower response 
K = -0.0274, 

= -70.88 
(section V.F) 


-0.0743 109.411 -2.32 0.0431 


with cross- 
coupling 
slower response 


-6.0791 199.841 -0.8662 1.0541 



and increased gain 
in tne coordinated 
branch = i.a, 

K = -0.2S7 
K = -70.88 
(section V.F) 

with cross- 
coupling slower 



response and 
modified pitch 
control law as 
section V.G 
= -0.0274 


-2.264 38.386 -0.00788 -0.6021 



Note: K^, K, K^p as referred in section V.F and V.G. 
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Figures 5.27 and 5.28 show the minimized inertial and 

kinematic cross-coupling, with feedback the 

angle-of-attack, with K. = -0.0274 and K = -70.88. 

1 

Comparing the values of the above minimized inertial and 
kinematic coupling, with the values of the kinematic and 
inertial coupling when the pitch channel has feedback the 
pitch angular rate (q) and = -0.0274, as given in table 
VI, it is concluded that the couplings when the pitch 
channel has feedback angle-of-attack (a) and the 
aforementioned particular compensator (Fgure 5.26), were 
reduced 228.29% for inertial and 190% for kinematic. 

Also, comparing the Figure 5.25 which illustrates the 
achieved inertial acceleration for circular airframe 
with minimized coupling effects and Figure 5.29 illustrated 
n„ for circular airframe with minimized coupling, it is 

LA 

concluded that as the coupling reduces, the coupling 
effects increases. 

H. CONCLUSIONS AND RECOMMENDATIONS 

The conclusions that follow are based on a single 
representative flight condition (Mach 3.95 and 60 kft 
altitude) and for the particular analysis which was done in 
this chapter. 

1. From the studies in this chapter, it is concluded: 
a. Transients in maneuver plane acceleration are 
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Fig. 5.27 



Minimized inertial Coupling (P.r) vs Time ( t) 
CBTT of Circular airframe 
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Fig. 5.28 



Minimized Kinematic Coupling (-P.6) vs Time ( t) 
CBTT of Circular airframe 
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R 




Fig. 5.29 Achieved Inertial acceleration (n^) vs Time (t) 
Modified pitch channel for minimization of 
coupling; CBTT of Circular airframe 
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.■v'» 



c 






caused by kinematic and inertial coupling between pitch and 
yaw dynamics through missile roll rate. Transients are in 
the form of overshoots and undershoots (section V.B) . 

b. Using as feedback in the nonlinear pitch 
channel for both airframes the pitch angular rate (q) , the 
transients in the maneuver plane exceeds the denoted 
requirements and the achieved maneuver acceleration is not 
acceptable (sections V.D and V.F). 

c. Using as feedback in the nonlinear pitch 
channel for both airframe, the angle-of-attack (a) and 
modifying suitably the pitch control law (Figure 5.5 and 
Figure 5.20), the transients are minimized below the 
desired requirement and the achieved inertial acceleration 
(n 2 ) is acceptable. (sections V.D and V.F). 

d. Using as feedback the rate of angle of attack 
(a) , the whole system becomes uncontrolable because the 
pitch tail incidence (5^) exceeds the limits. (sections 
V.D, V.F) . 

e. Using as feedback the rate of angle of attack 
(a) and modifying suitably the pitch control law (Figure 
5.26), the kinematic and inertial coupling are reduced 
dramatically, comparing with their values when the used 
feedback is the pitch angular rate, (sections V.D, and V.G) 

f. Using as feedback the rate of angle of attack 
(a), in the modified pitch channel (Figure 5.26) for 
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minimization of kinematic and inertial coupling, the whole 
system becomes uncontrolable , because the pitch tail 
incidence ( 6p) exceeds the limits, (sections V.E and V.G) 

g. Minimizing the inertial and kinematic 
cross-coupling effects, the values of kinematic (-pS) and 
inertial (pr) coupling increase, (section V.E and V.G) 

h. The above conclusions show, that for the 
particular airframe, flight conditions and analysis is not 
possible to achieve minimization of the kinematic and 
inertial cross-coupling effects and simultaneously 
minimization of the cross-couplings themselves. 

2. Recommendations 

Considering the studies of this work and the 
conclusions, it is recommended: 

a. Further analysis of the minimization of the 
kinematic and inertial cross-coupling effects and the 
cross-coupling themselves in other airframes, and with 
other analysis to achieve simultaneous minimization. 

b. Further analysis of the CBTT autopilots in 
other flight conditions for verification of the above 
results. 



246 



I 



APPENDIX A 



MISSILE SIZING AND MASS PROPERTIES 
In order to provide a realistic missile based on the 
configuration concepts tested aerodynamically in [Ref. 1] , 
the models have the following characteristics. 

Characteristic Circular Elliptical 



Length (in) 


168 


168 


Max. Diameter 


24 




Max. Major axis (in) 




41.57 


Max, Minor axis (in) 




13.86 


Weight (lbs) 


2525 


2475 


Ixx (slug-ft^) 


40 


110 


lyy (slug-ft^) 


804 


790 


1^2 (slug-ft^) 


810 


853 


c.g distance from L.E (in) 


100.8 


100.8 


Reference length d (ft) 


2 


2 


2 

Reference area S(ft ) 


T7 


IT 
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APPENDIX B 



NONLINEAR AERODYNAMIC DATA 

The aerodynamic data were taken from [Ref. 2] . For 
reference, the normal and pitching moment curves at M = 
3.95 have been produced in fig. B.l and B.2 for the two 
configurations in [Ref. 2]. The aerodynamic derivativies 
C , C , with respect to sideslip angle b, yaw control 
& , and roll control are presented in figures B.3, B.4 
and B.5 of [Ref. 2] as used in the computer simulation, 
namely, as piece wise linear segments for ease in 
interpolation . 
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J^PPENDIX C 



REQUIREMENTS FOR UNCOUPLED AUTOPILOTS 



The requirements for the classical design technique of 
the uncoupled autopilots are given in [Ref. 2] and are the 
following ; 

1. Acceleration time response 

a. 63 percent time constant of 0.5 seconds for a 
step command of acceleration at the flight condition of 
interest, (i.e M = 3.95 and altitude 60kft) and small 
angles-of-attack. This response is representative of a 

tactical missile of this size. 

b. Overshoot 10 percent 

c. Zero steady state error in acceleration to 
reduce variations of guidance navigation gain. 

2. Relative stability 

Gain margins 6db, phase margins 30 deg with a 
goal of 12 db and 50 deg. 

3. High F requency Attenuation in Actuator Command 
Branch 

It must be ^ 15 db at 100 rad/sec and zero 
angle-of-attack. This requirement will provide sufficient 

high frequency attenuation for 30 Hz actuator and for body 
bending modes, when high frequency filters are added. This 

requirement also limits autopilot speeds of response. 
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CCMMANO FOR CIRCULAR AIRFRAME 
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